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1.0 INTRODUCTION 


The Large Antenna Multi-channel Microwave Radiometer (LAMMR) is a 
mechanically scanned multi -frequency passive microwave sensor that re- 
ceives power emanating from the surface of the earth in seven spectral 
channels for both horizontal and vertical polarizations. The seven 
microwave frequency bands lie between 4.3 and 36.5 GHz and are extracted 
using a single high-gain antenna which scans the earth at one revolution 
per second, mapping its surface and atmosphere. The geophysical proper- 
ties of sea surface temperature, oceanic wind speed, precipitation, water 
vapor, sea ice concentration, snow and land presence can be inferred from 
the antenna brightness temperatures. The finite parabolic antenna 
receives significant energy within .26 to 1.3 degrees depending upon 
microwave frequency, but also receives some non-negl igible side lobe 
energy . The LAMMR embodies a 4.0 meter aperture antenna which scans the 
earth in a circle of 43.6° half cone angle from nadir. Data sampling 
occurs during the forward 150° of rotation as shown In Figure 1.1. Table 
1.1 shows the major LAMMR characteristics. 

LAMMR software specifications have been written for LAMMR ground 
processing. There is a need to determine more computationally-efflcient 
antenna temperature correction methods in compensating side lobe contri- 
butions especially near continents, islands, and weather fronts. In a 
previous study/ antenna patterns were shown to be digitally modifiable 
and alternative philosophies were established for Antenna Pattern 
Corrections (APC). One of the major realizations of this study was that 
APC processes did not accomplish the implied goals of compensating for 
the antenna side lobe influences on brightness temperature. A-priori 
knowledge of land/water locations was shown to be needed and had to be 
incorporated in a context sensitive APC process if the artifacts caused 
by land presence is to be avoided. The high temperatures in land regions 


"•■Alternative Formulations of Antenna Pattern Correction For 
N0SS/LAW1R Real Time Processing, BTS Report 8TS-FR-81-153 under contract 
No. NAS 5-26280. 
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PERFORMANCE 


• Sea Surface Temperature accurate to within - 1.5* at 25 km resolution 

• Wind Speed 0-50 m/s, accurate to within - 2 m/s or 10?, whichever is greater 
at 17 km resolution 

• Sea Ice Concentration to within ± 15? at 9 km resolution 

- Type classification - first-year, multi-year 

- Age - inferred thickness to within * 2 m 

• Atmospheric Water Vapor to within * 0.2 g/cm 2 at 9 km resolution 


LAMMR TECHNICAL CHARACTERISTICS 


• Orbital Altitude 

• Weight 

• Digitization 

• Data Rate 

• Power (nominal -average) 

• Frequency (GHz) 

• Wavelength (cm) 

• Antenna Aperture 

• Polarization 

• Cross-Polarization Isolation 

• 3 d8 Beamwidth (deg) 

• Beam Efficiency 

• Sensitivity Ta (K) 

• Calibration Accuracy (K) 

• Scan Nadir Angle 

• Viewing Zenith Incidence Angle 

• 2 -Axis Viewing Angle Stability 

• Footprint (km x km) 

• Scan Rate 

• Active Scan ARC (Forward) 

• Swath Width (150* Scan) 

• IFOV's Sampled Per Scar 

(Each Channel) 


700 km 
320 kg 
12 bits 
54 kb/s 

350 w (Nominal Average) 

4.3 5.1 5.5 10.65 18.7 21.3 36.5 

7.0 5.9 4.5 2.8 1.6 1.4 0.8 


- - 4.0 m ii. — 

HORIZONTAL ANO VERTICAL 
— 17 db minimun 


1.3 

1.1 

0.9 

0.6 0.3 0.25 

0.25 






0.6 

0.7 

0.8 

0.8 0.8 0.3 

1.5 

1.0 

1.0 

1.0 

1.0 1.5 1.5 

1.5 


43.5° 

50.0* 

0 . 02 * 

23x36 20x31 16x25 11x17 5x8 5x7 5x7 

1 Hz 

150* 

_____ — ___ 1360 km 

256 } { 512 


Table 1.1 LAMMR Baseline Technical Characteristics 
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can severly bias the lower ocean response even though the side lobe gain 
through which land and weather fronts are viewed may only be a few 
percent of the main beam gain. Linear context-free techniques fall short 
of achieving this goal. It is for that reason that this study of antenna 
pattern corrections and side lobe compensation was Initiated to extend 
the results of the previous effort. 

1.1 APC Processing Goals 


There are two major goals associated with the context sensitive 
formulation of antenna pattern correction and side lobe compensation: 

• An error analysis and assessment to show that the proposed 
context sensitive schemes are satisfactory, and 

• Demonstration that acceptable error allowances can be 
achieved within reasonable computational loads. 

In the error assessment area, a least squares error measure was defined 
under some reasonable underlying statistical assumptions concerning the 
signal characteristics relative to noise. The use of a priori knowledge 
in the processing of antenna temperatures compounds the issue with regard 
to error measures which this study resolves. A unified scheme and/or 
philosophy of processing was established. We highly recommend a context 
sensitive APC process which we believe yields a good tradeoff between 
contending factors and which significantly reduces the computational 
burden. This report demonstrates its effectiveness. In the 
demonstration area, the computation of APC coefficients under alternative 
signal -to-nolse conditions, under various matrix truncation conditions, 
and under aliasing and non-aliasing sampling spaclngs wore selected and 
applied to the unmodeled portion of the antenna temperature field. All 
analyses and supporting rationale were completed and the demonstration of 
their effectiveness using a synthetic Image model was accomplished. Work 
was initiated on trying to apply the technique to SMMR data, but time did 
not permit us to complete this extention. 
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1.2 Error Assessment 


Within our previous study an error assessment was made Incorporating 
many error measures: 

• Error variations due to geometry within the swath, 

• Noise amplification resulting from the APC process, 

• Alternative error measures associated with different signal 
auto-correlations, 

• APC matrix truncation error, 

• Error in being able to achieve correct brightness tempera- 
tures as a result of sensor noise, or data sampling limi- 
tations (aliasing errors), and 

• Error reduction due tc Incorporating a priori knowledge In 
a context-sensitive APC process. 

Least squares error estimates were made under all of the alternative 
conditions, but little use can be made of these estimates in a context- 
free processing scheme since environmental conditions vary greatly in the 
real world. 

Context sensitive APC processing adapts to varying environmental 
conditions. It Incorporate a-priori knowledge about land mass, ice, and 
water locations and models the brightness temperatures In those areas. 
This brightness temperature model Is used as a bareline from which an 
antenna temperature can be computed and compared with measured antenna 
temperatures T ft . The variation of T^ from the model can be 
corrected to variations of brightness temperature estimates <5Tg from 
the brightness temperature model so that a final estimate of brightness 
temperatures can be made. Also, varying statistical correlations can be 
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used to reflect coherences that may be different over land than they are 
over water, etc. The main advantages of this approach are: 

• Gibbs artifacts generated by correcting the finite spectral 
response of the antenna near land/water boundaries are 
overcome. 

• The computational system Is consistent In that no special 
attention has to be given to those pixels near each end of 
the swath traces as would otherwise be required and which 
is described in the algorithm freeze report. 

• It utilizes more realistic statistical assumptions in 
establishing a matched filter for modulation transfer 
function compensation which can be different over different 
fields. 

• The algorithm is computationally more efficient than the 
context-free cases. 

In the theory, we found no need to deviate from a least squares 
error assessment model except to realize that the underlying statistical 
assumptions may not be context sensitive; differing over land, water, 
snow and ice. 

1.3 Real Time Processing Constraints 

Within the BTS study, concepts for reducing computational workloads 
are suggested. Some of these concepts are: 

• Reduction of matrix size while meeting acceptable error 
allowances for modulation transfer function compensation of 
antenna temperature variation from the model. 
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• Use of world surface map features for adaptively processing 
the data, 

• The use of a priori brightness temperature models for data 
Interpretation, and 

• Application of APC process In Instrument coordinates. 

This final report clearly identifies factors and assumptions used in 
the proposed schemes. It is recommended that APC be performed In instru- 
ment coordinates so that no geographically-referenced binning be required 
as in SMMR. Computer configuration requirements such as core, data 
representations and processing speeds are presumed, other factors such 
as programming language, processing structure, algorithm descriptions, 
mathematical specifications, and data structures for both input and out- 
put are consistent with the LAMMR algorithm freeze specifications. 

Our report Is organized as follows. Section 2 contains our execu- 
tive summary. Therein contained is a description of the context 
sensitive APC philosophy, the presentation of computational results, a 
discussion of pertinent findings, and a few recommendations. Section 3 
consists of the analysis conducted during the course of this work. It 
also includes examples of some of the compu, '■ions which were conducted. 
Related analysis needed for understanding this section are to be found in 
an appendix. Section 4 describes our software code In modular form and 
references several appendices containing the detailed runs and code 
listings. Section 5 gives our conclusions and recommendations, it 
identifies problem areas and delineates approaches to their solutions. 
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2.0 EXECUTIVE SUMMARY 


This section describes the essential results of the context 
sensitive Antenna Pattern Correction (APC) and Side Lobe Compensation 
(SLC) study. It contains the objectives the rationale for selected 
philosophies, study approach, results and recommendations. 

2.1 APC and SLC Objectives 

There are two major objectives In this study. 

• Alternative approaches to context sensitive antenna pattern 
corrections and side lobe compensation that are considered 
to be computationally viable are to be formulated utilizing 
a statistically rigorous derivation and geographic model 
information so as to establish alternative philosophies 
and/or assumptions needed in these theories. This analysis 
is to be an extension of work that had been previously 
performed. 

• To perform some computational demonstration that are 
sufficient to establish the algorithmic feasibility of 
context sensitive antenna pattern corrections and to make 
recommendations for an approach to be used In LAMMR 
processing. 

2.2 Rationale for Selected Philosophies 

In previous work by BTS and others, the words "antenna pattern 
correction" implied a process by which antenna temperatures are converted 
to estimates of brightness 'emperatures over some spatial field. Linear 
convolusion processes were the only processes previously considered for 
this conversion. In our previous study. It was shown that the coeffi- 
cient to be used in this estimate of brightness temperature can be 
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rigorously derived as long as there Is knowledge of the brightness 
temperature auto-correl ation function. That is, one requires statistical 
knowledge of the fluctuations to be expected In the brightness tempera- 
ture field. In that study It was shown that the presence of land in the 
brightness temperature field causes a disturbance which Is not adequately 
modeled by an^ reasonable statistical auto-correlation function 
assumption. Therefore, the underlying theory of linear convolutions for 
brightness temperature estimates are not suitable in the land/water 
boundary regions. Since It is precisely the coastal regions where 
accurate knowledge of brlgntness temperatures have their greatest 
utility, it becomes essential to solve the dilemma. In so doing, we have 
found it convenient to separate the brightness temperature estimation 
process Into parts. The first part we call side lobe corrections , and 
the second part we call antenna pattern corrections. 

It was previously thought that antenna pattern corrections would 
solve the problems associated with side lobe influences. It was known 
that when a bright Island appeared under one of the major side lobe 
peaks, that brightness temperatures would radiate into the antenna there- 
by affecting •* ts antenna temperature measurement; an erroneous reading in 
the surrounding waters is created. Since antenna pattern corrections 
modify the side lobe structures of the effective antenna pattern, it was 
thought that these influences could be subdued numerically and indeed 
such Is the case, except that this reduction causes a significant widen- 
ing of the main beam with a subsequent loss of resolution, i.e., the 
process of reducing side lobe influences blurs the Imagery. Antenna 
characteristics already cause a blurring of the imagery; rapid spatial 
variations or high spatial frequencies are severely modulated as the 
antenna radiometer responds to these variations. The typical engineering 
process of modulation transfer function compensation would digitally 
amplify the higher frequency components, so as to compensate for their 
loss as a result of the antenna system characteristics. It is easy to 
show, however, that such a compensation which amplifies higher spatial 
freqencies will narrow the main beam of the effective antenna pattern but 
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simul taneously introduce a large number of relatively high amplitude 
(positive and negative) side lobes. The existence of a small island 
and/or a land water boundaries creates significant artifacts in the 
brightness field estimated by such an effective antenna pattern. 
Approaches which are based on optimum estimates of brightness temperature 
fields indicates that the modulation of the spatial response of antenna 
system should be compensated as previously described in a well known 
manner known as Weiner Filtering. However, this analysis must utilize an 
a priori estimate of the autocorrelation in brightness temperature 
fluctuations about some expected mean brightness temperature. The usual 
assumption that these brightness temperature fluctuations are spatially 
independent is valid in uniform regions, such as, over land or over 
ocean. That assumption is totally inadequate in the vacinity of 
land/water boundaries and the correction process generates artifacts. To 
avoid the artifacts, land presence must be explicitly or implicitly 
recognized. Several approaches were considered for incorporating our j* 
priori knowledge of land/water differences in side lobe compensations. 

• The geometric knowledge of land/water interfaces could be 
utilized to generate an a priori field of expected 
brightness temperatures over the surface of the Earth. 

• The same model could be parameterized with two brightness 
temperatures; a parameter for the brightness temperature 
over water and a separate parameter for the brightness over 
land. 

• An implicit approach could be taken in which areas above a 
set threshold would be considered land while those below 
the threshold would be considered water. In the implicit 
scheme land /water boundaries are extracted directly from 
the antenna temperature measurements. 
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As a result of ou' analysis it was concluded that the second of these two 
schemes was appropriate for LAMMR. The reasons for this choice are as 
f ol 1 ows : 

An unparameterized model as in the first bullet would be expected to 
have erroneous a priori brightness temperature values over land and water 
within the model. The differences between these assumptions and the 
actual values would lead to a spatial discontinuity at the boundary 
causing some residual artifacts. The third bullet was rejected because, 
in the processing algorithm for LAMMR, knowledge of land/water locations 
are required outside the swath of observed data. A simple thresholding 
technique would not yield this information and so a world boundary map 
would be required. This requirement would cause the algorithm for pixels 
along the boundary of the swath to be different than those on the 
interior of the swathe complicating, therefore, the SLC process. Hence, 
it was decided that the second approach of a parameterized model 
utilizing a priori knowledge of land/water locations was most appropriate 
for LAMMR. 

We must now show how such a model can be used to remove the influ- 
ences of the side lobes; the process is called side lobe compensation 
(SLC). The side lobe compensation process is recognized as one which has 
an influence only in the presence of land; large oscillating side lobes 
with positive and/or negative contributions to the antenna temperature do 
not generate artifacts when the field of observation is relatively 
uniform. Side lobe influences can be removed in two ways. 

• Removal of the side lobe structure and antenna pattern 

• Removal of the influences caused by land 

As previously described, the reduction of side lobes by affecting 
the antenna itself has two major disadvantages. First it causes the main 
beam which is already wide to be further widened, thereby, leading to 
image blur and loss of resolution. Secondly, the process of modifying 
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the antenna pattern to remove side lobes lead to large convolution 
matrices which are computationally expensive and needlessly complicates 
the system. Through modelling It Is possible to remove the Influences 
that the land/water boundaries have on brightness temperature estimates. 
In this process the parameterized model of brightness temperatures Is 
convolved (a one time process) with the antenna pattern of the instrument 
generating a model antenna temperature field. This antenna temperature 
field is compared with the observed data to determine the values of the 
to brightness temperature parameters which best explain the observed 
data. Now differences between the observed antenna temperature and the 
modelled antenna temperatures are antenna temperature unexplained 
variation. These unexplained variations are no longer influenced by the 
presence of land and/or water; the gross mismatched between brightness 
temperatures over land and water have been removed. By moving to the 
space of unmodelled antenna temperatures, we have eliminated the 
influences of land beaming into he side lobes of the antdnna pattern and 
have therefore performed a side lobe compensation function. Hence, we 
call the process of moving to the space of unexplained antenna tempera- 
ture variations the process of side lobe compensation. Note that we have 
not yet corrected for the response function of the antenna system which 
is appropriately called antenna pattern corrections . The antenna pattern 
correction process is applied to the unexplained antenna temperature 
field to generate estimates of an unexplained brightness temperature 
field. This process must simply compensate for the loss of high 
frequencies that have been significantly modulated by the observation 
process. In our previous report, the methodology for linear APC 
processes have been e v jstively discussed. It was shown that most of 
the high frequency compensation process can be accomplished with a 
relatively small matrix. A 5*5 array of coefficients is adequate and it 
is even possible that a 3*3 APC process would be just as good. It must 
be recognized that the set of coefficients that are used for high 
frequency compensations do depend upon a priori assumptions in the 
spatial coherence of the urnnodelled brightness temperature field. As a 
first step in this assumption, the unmodelled brightness temperature 


13 



Business and Tecmnoukucal Systems, Inc 


field would be assumed spatially uncorrelated. This assumption can be 
tested over land, water, ice and snow once real data is obtained. 
Modifications to APC matrix coefficients could then be modified context- 
urally and applied to these separate categories. This refinement, 
however, is considered to be unessential as an initial system 
specification. 

2.3 Study Approach 

The approach used to accomplish the task associated with the 
context-sensitive antenna pattern correction and side lobe compensation 
study are embodied in Figure 2.1. Here the functional milestone chart of 
the antenna pattern correction study is presented for each of five tasks 
and associated subtasks. Also shown are tasks which have been completed 
in a previous study for context-free APC evaluation. 

The tasks which were herein completed include: 

• Analysis and planning 

• APC benchmark software development 

• APC weighting coefficients 

• Verification testing and analysis 

• Documentation and demonstration 

In the analysis and planning tasks the APC/Data Bank Algorithms 
which are described in the LAMfitf freeze report dated October 1980 were 
analyzed and checked for completeness, accuracy, adequacy and 
cumbersomeness. From our previous work it was determined that the 
algorithm as described had elements of a context-sensitive nature which 
had to be applied along the boundaries of the swath. In particular, 
land/water flags were needed in the region just outside the swath of 
available data so as to be able to perform an antenna pattern correction 
function for several pixels near the beginning and end of each scan 
line. This need for land/water flags coupled with the realization that 
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Figure 2.1 Functional Milestone Chart of Antenna Pattern Correction Study 
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the matrix multiplication process may not compensate for the influences 
of side lobe interferences has lead us to consider a unifying context 
sensitive APC process for LAMMR. The use of the land/water flags are 
needed to interpret the observed antenna temperatures so that artifacts 
free brightness temperatures can be drived from them. A unifying 
philosophy was established which treated the beginning and end of scan 
data no differently than the rest of the data and which took advantage 
everywhere of our knowledge of land/water boundaries. Software was then 
planned for developing a simulated scene and trade-off evaluation plans 
were established for ut^Ting this scene to demonstrate this context 
sensitive approach. A work study plan was generated. 

In the APC weighting coefficients tasks, context-sensitive 
algorithms were studied and analyzed. The work un the generation of 
coefficients for antenna pattern corrections was extended so that we 
could better understand the trade-off between the modulation transfer 
function compensation and desires for reducing the amount of arithmetic 
needed in the performance of this function. Studies were conducted 
utilizing alternative sets of matrix coefficients for modifying the 
antenna pattern and the differences were evaluated in the context of the 
LAWiR Instrument. 

In the APC benchmark software development task, software was 
developed to synthetically simulate a scene consisting of land and 
water. This scene consisted of a rectangular island completely 
surrounded by water. The geometry of this island is such that the 
effects of the context-sensitive algorithms could be studied both in the 
vicinity of straight boundaries as well as along sharp curves. Software 
was also generated utilizing the synthetic scene to simulate the effects 
of the antenna viewing this scene. This simulator generated antenna 
temperatures from the given brightness temperatures of the land and 
water. 


16 



BUSINESS AND TECHNOLOGKAL SYSTEMS, INC 


In the verification testing and analysis task, both context-free and 
context-sensitive antenna pattern correction and/or side lobe 
compensation algorithms were applied to the simulated scene. The results 
verified that context-sensitive algorithms outperformed context-free ones 
especially In the vicinity of land/water boundaries. Plans were made to 
apply these algorithms to SMMR data. A world scene was extracted for use 
as the land/water flags to replace the synthetic scene. Studies were 
made of available SMMR data which could be applied in a complex 
water/land environment. It was decided that the region of the 
medlteranian which Included Italy would be Ideally suited to show the 
advantages of context-sensitive processing. However, several problems 
arose that did not permit us to complete the simulation on SMMR data. 
These problems Include: 

• The Government furnished land boundary flaged software 
which was was coded in such a manner that the region 
between water and land had mixed pixels designations. This 
would have been okay if such a mixed pixel was a single 
pixel. However, for some reason this file blurred the 
land/water boundaries by as much as ten pixels so that the 
region of interest lacked the resolution needed to demon- 
strate the utility of our algorithms. 

• In an effort to obtain raw antenna temperature measurements 
as observed by the SMMR instrument it was ascertained that 
such raw data did not exist in a format that would be easy 
enough for us to use. The only data set that was conven- 
iently available for processing had been already "binned" 
and had therefore, totally destroyed relationships from 
which an antenna pattern correction function could be 
expected to give compensated results. 
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Efforts, however were extended to remedy these situations. For 
example, the land/water boundary flag algorithms were modified in an 
attempt to avoid the bluring function that is embeded in that algorithm. 
This proved to be a difficult task because there is no documentation of 
the raw data from which the land boundary information is extracted nor is 
there any documentation within the code which define the variables and 
the functions which are performed therein; reading the code makes one 
wonder what is going on and why. Efforts to modify this code to generate 
sharper boundaries have thus far failed and is perhaps the reason why the 
boundaries are blurred in the first place. Some effort was extended to 
understand the nature of available SMMR data. It was ascertained that 
the binning operation completely swamped the effects of the SMMR antenna 
pattern so that it was dubious that the SMMR context-sensitive simulation 
would prove anything useful for LAMMR. Because of difficulties that 
would be encountered in a proper simulation this effort was abandoned in 
favor of a more extended study of the simulated scene. As a result of 
work done on the simulated scene some trade-off evaluations were 
conducted which lead to the algorithm proposed in this documentation. 

In the documentation and demonstration subtask context-free and 
context-sensitive software was documented. Final report was prepared and 
demonstration material was generated as required. This task also 
included the generation of monthly progress reports. 

2.4 Study Results 

A compilation of the essential results of our study follows: 

* Alternative Approaches : Detailed analyses of the context 

sensitive antenna pattern correction and side lobe compen- 
sation processes indicated that some model assumptions and 
some statistical assumptions are absolutely required. 

Different assumptions in this detailed analysis 
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leads to alternative philosophical approaches to antenna 
pattern correction as well as to alternative approaches to 
the side lobe compensation functions. 

• Alternative APC Coefficient Computations : For each of the 

alternative philosophies, additional approximations can be 
made in the computation of antenna pattern correction sets 
of coefficients. A matrix approach has been developed when 
the set of correction coefficients are truncated to a 
finite set. Fourier techniques are applicable when the 
antenna pattern rotation effects can be locally Ignored. 
Techniques for computing coefflents using the Fourier 
approach have been demonstrated for both over-sampling and 
under-sampling conditions, and have been shown to agree 
quantitatively with the matrix truncation approach. In the 
context-sensitive approach, it was ascertained that the APC 
function is one of compensating for the modulation transfer 
function of the system and that this can adequately be done 
using a 5 x 5 matrix approach. 

• Instrument Coordinates : It has been shown to be computa- 

tionally efficient to perform the computations in instru- 
ment coordinates in that the raw data exists in this coor- 
dinate sysem. Resampling of the raw data to fit a priori 
selected geographic blocks on the surface of the earth as 
in SMMR processing hae been found to be unrealistic and 
undesireable. Hence, "Binning" should be avoided. 

• Data Interpretive Techniques : A data interpretive tech- 

nique has been developed whereby the discrete set of 
antenna pattern correction coefficients can be computed by 
solving a truncated system of linear equations. Because 
the antenna pattern varies significantly over the distance 
associated with a sample, it is necessary to retain a con- 
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tinjum representation of antenna patterns for the computa- 
tion of matrix elements which enter Into the system of 
linear equations for APC. 

• Geometric Model : A complete geometric model which relates 

the Instrument In nominal spacecraft coordinates to a 
geographical ly-referenced coordinate system with associated 
scanning geometry has been developed. Geographic 
referencing is needed to extract the land/water flags. 

• K Space : Because of the finite aperture of the LAMMR 

antenna, it is strictly correct that the Fourier transform 
of the antenna pattern has a sharp spatial frequency cut- 
off. It has been determined that this statement is correct 
only in the K Space corresponding to the projection of the 
Poynting vector onto a plane orthogonal to the antenna 
axis. Spatial transforms from instrument coordinates to K 
Space have been derived so that Fourier transforms in K 
Space can be utilized, taking advantage of the share fre- 
quency cut-off that occurs in that space. This is useful 
for generating consistent antenna patterns and for 
performing convolutions with simulated scenes. 

• Earth Rotation ; It has been determined that earth rotation 
effects are significant in two respects. First, the ground 
velocity vector and the instrument axis corresponding to 
the center of the sampling angles do not coincide as a 
result of earth rotation; and second, the separation 
between lines alters as the ground velocity is modulated by 
earth rotation effects. 

• Antenna Pattern Representation ; It has been determined 
that adequate representation of antenna patterns are viable 
utilizing linear interpolation if four points per nominal 
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sampling interval are used to represent the antenna pattern 
with sequences of straight line segments. All simulated 
results utilize this approximation when convolution 
techniques are used. 

• Error Measures : Several error measures have been defined 

in the context of this study. However, the least squares 
difference between the estimated brightness temperature and 
the corrected brightness temperature has been selected as 
the only reasonable criteria for measuring the effective- 
ness of alternative antenna pattern correction processes. 

It is realized, however, that this same error measure can 
be applied to alternative models. A discussion of an 
appropriate mix of models is to be found in our previous 
report. 

• Noise Amplification : Each of the antenna correction pro- 

cesses enhances frequencies (or suppresses frequencies) 
that were modulated by the antenna pattern of the LAMMR 
instrument. As a result of the spatial frequency enhance- 
ment, noise at those frequencies are likewise enhanced so 
that a total merit of noise amplification is pertinent to 
the analysis. Codes for computing noise amplification have 
been developed and used in generating optimum ARC matrix 
coefficients. 

• Effective Power Spectra : In the truncated matrix approach 

to antenna pattern corrections, a linear combination of 
local data is chosen as a means to enhance the spatial fre- 
quencies of the unexplained antenna temperatures that have 
been modulated by the antenna. As a result of this linear 
combination, the resulting brightness temperature will have 
a corresponding point spread function and associated power 
spectra. Codes for computing both of these curves in both 
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one- and two-dimensions have been developed. The large 
side-lobe oscillations In these patterns are not Important 
when applied to the "unexplained" antenna temperatures 
because the effects of land have been removed and can no 
longer have a large Influence on the brightness temperature 
measurement. 

• Aliasing Versus Signal -to-Noise : From our analysis, we 

find two competing sources of error which govern the choice 
of antenna pattern correction coefficients. These compet- 
ing factors are signal -to-noise ratio and the degree to 
which the data has been unoersampled. It has been demon- 
strated that when undersampling is such that the signal 
strength at the Nyquist frequency is larger or comparable 
to the noise power, then the affects of aliasing and errors 
associated thereto completely dominate the antenna pattern 
correction process independent of the signal -to-noise 
ratio. When the context sensitive SLC processes are 
invoked, the effects of aliasing do not appear to be 
significant. 

• Truncation : In order to keep the number of computations in 

the antenna pattern correction process reasonable, it is 
necessary to truncate the array of coefficients to he 
applied to the data in converting antenna temperatures to 
brightness temperatures. It has been shown that there is 
little to be gained by considering large arrays for antenna 
pattern corrections when SLC processes are used. 5x5 
sets of coefficients are recommended, but 3x3 
compensation appears adequate. 

• Undersampl ing : A model of the antenna pattern associated 

with a four-meter dish, its electronics, its motion and 
sampling strategy, as planned for the LAMMR instrument has 
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been developed. Although the three db points agree In gen- 
eral with those previously published for this Instrument, 

It has nonetheless been determined thit In the context of 
reasonable signal -to-nolse ratios, and In the context of 
our ability to enhance spatial frequencies that had been 
modulated significantly, that the planned sampling strategy 
significantly undersamples the data over most of the fre- 
quency channels. As a result, aliasing Is expected to 
occur, and limit the spatial resolution which could other- 
wise be acquired. Techniques for dealing with aliasing 
have been developed; aliasing Is dominant In most cases but 
may not be significant when geometric models are used. 

• Gibbs Phenomena : The discontinuity which occurs at the 

land/water Interface contains unbounded spatial frequen- 
cies. Because of the finite dimensions of the antenna, 
only those spatial frequencies that are less than the cut- 
off frequency will be sensed by the antenna. Under Ideal 
conditions where noise Is Insignificant relative to all or 
most of the observed spatial frequencies, compensation for 
the modulation of each frequency can be digitally Imple- 
mented. Even If this were done, however, the resulting 
edge woulo show the characteristic "ringing" associated 
with Gibbs phenomena. It has been determined that addi- 
tional a priori knowledge associated with the natu e of 
land/water discontinuities must be utilized and Incorpor- 
ated within the antenna pattern correction process In order 
to avoid the Gibbs ringing. A method similar to one util- 
ized by Dr. Schell for radio astronomy has been developed 
for this purpose, and is designated as a context-sensitive 
antenna pattern correction process. This process In effect 
compensates for the side lobe Influences of land and is 
called side-lobe compensating (SLC). 
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• Context-Sensitive Processes : Context-sensitive antenna 

pattern correction processes have evolved during the course 
of these studies. One very simple process is to utilize a 
data base consisting of the difference between a modeled 
antenna temperature and Its associated brightness tempera- 
ture as a means of correcting observed antenna temperatures 
which appropriately incorporates the a priori knowledge of 
the location of land/water boundaries. A parameterized 
version is the one recommended for LAMMR. 

• World Brightness Map : As a result of our investigation, it 

appeared feasible to utilize a world brightness map con- 
taining a priori knowledge of expected brightness tempera- 
tures over the surface of the earth. A model of the 
antenna could be applied to this world brightness data 
base, generating expected antenna temperatures associated 
with this model. Differences between antenna temperatures 
computed In this way and the modeled or statistical ly- 
observed brightness temperature constitutes a bias which 
should be utilized for antenna pattern correction. This 
bias has been determined to consist of two parts; first, a 
part resulting near discontinuities in brightness tempera- 
tures such as land/water interfaces, and second, biases 
resulting from the tails of the antenna pattern resulting 
from the amount of land to be found in the far tails of the 
distribution. These biases may have been observed in SMMR 
data, and may have beer, erroneously attributed to radio 
frequency Interference. Attempts are being made to 
generate a World Brightness Map suitable for LAMMR 
processing. 
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2.5 Recommendations 


There are two major recommendations to be made resulting from this 
study on alternative processes of context-sensitive antenna pattern 
correction and side lobe compensation. These recommendations are: 

• That a parameterized context-sensitive antenna pattern 
correction process be implemented for LAMMR in that no 
single context-free process will adequately compensate for 
the antenna modulation of higher spatial frequencies while 
simultaneously avoiding the artifacts and errors generated 
by bright spots, glaring into the side lobes of the antenna 
pattern. 

• That the LA>WR APC algorithm be separated into two 
functions. First is a side-lobe compensation function in 
which a parameterized World Brightness Hap is used to 
estimate the expected t^nna temperatures and an 
"unexplained" antenna temperature field is generated as a 
representation of the data. Second, a matrix method of 
compensation be applied to correct for the modulation 
transfer function of the antenna system. 
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3.0 ANALYTIC RESULTS ^CEDING PAGE ELA .\'K NOT FUMED 

Analysis has been performed on the context sensitive antenna pattern 
correction and side lobe compensation problem. In the previous study two 
characteristically different types of antenna correction processes have 
been identified; context-free and context-sensitive antenna pattern 
correction processes. Very careful analyses have been conducted so as to 
characterize the problem in the continuum as opposed to the discrete 
domain. In order to do this, BTS has developed a data interpretive 
technique in which exact continuum representations are used in 
conjunction with integral representations of cost functions, in order to 
establish optimum coefficients in the matrix estimation of brightness 
temperatures as seen through some desired "eye“ window function. The 
technique is a generalization over that used by A. Stogryn 1 in 
establishing brightness temperatures from scanning radiometer data. The 
procedure involves many steps as indicated by the detailed analysis which 
is presented in Appendix A. 

In this section, we attempt to illustrate the concepts of context 
sensitive approaches to Brightness temperature estimates and leave any 
new mathematical analysis to Appendix A as much as practicable. 

3.1 Context-Free Versus Context-Sensitive Processes 


The distinction between a context-sensitive and a context-free 
process of antenna pattern correction is easily seen. All context-free 
antenna correction processes involve the computation of a linear 
combination of observed antenna temperatures in estimating the brightness 
temperature at any point. This 1 inear combination may be performed via 
matrix multiplication or other approaches such as the use of Fast Fourier 
Transform techniques. The important distinction is that the coefficients 
that are used in the linear combinations are independent of antenna 

*IEEE transaction on antennas and propagation Volume 18-26, No. 5, 
September 1978. 
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temperature measurements or patterns that are observed in the data, nor 
are the coefficients dependent upon other a priori knowledge such as the 
location of land/water boundaries. On the other hand, a context- 
sensitive process depends upon a priori knowledge and is conceptually 
more complex, although not necessarily computationally less efficient. 

Context-Free Theory 


The coefficients to be used in a linear context-free process are 
predetermined coefficients independent of observed brightness tempera- 
tures and independent of other knowledge of the real earth brightness 
temperatures. This does not imply that these coefficients are indepen- 
dent of any a priori knowledge. Indeed, it is necessary to utilize a set 
of statistical assumptions (an a priori model ) in deriving the desired 
set of coefficients to be applied in the APC process; different context- 
free correction philosophies will lead to different sets of coefficients. 
The following types of assumptions are needed for a context-free theory: 

• Choice of cost functions, 

- Least squares 

- Minimum variance 

- Minimum entropy 

• For the case of minimimi variance, a polynomial model of bright- 
ness temperatures may be local 1 required 

- Mean model 

- Constant radiant 

- Higher order splines, etc. 

• A priori knowledge of either the brightness temperature auto- 
correlation f notion or the auto-correlation function of the var- 
iance about a polynomial model is required 

- Zero correlation model 

- Exponential correlation model 

- Linear auto -correlation function model. 
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As can be seen from the above list of potential context-free models, 
many alternative philosophies are possible, and each philosophy will lead 

to a separate approach to antenna pattern corrections, yielding different 

sets of coefficients. In addition to the above assumptions, it is also 

possible to include other factors such as a penalty function of the 

Backus-Gilbert type which could be utilized to help suppress side lobes 
in exchange for minor degradation in the spatial resolution of the main 
beam. In the computer simulations that were performed, very specific 
philosophies for context-free and context-sensitive theories were 
chosen. The context free philosophy was oriented toward the 
establishment of maximum accuracy of inferred brightness temperatures 
values in open ocean. Here, a minimum variance approach was chosen in 
which the geometric model of brightness temperatures was assumed a priori 
to be a constant. Fluctuations in brightness temperature from this 
constant was assumed to have no spatial correlation whatsoever, but with 
a mean fluctuation signal of S . The a priori power spectra therefore 
was S 2 and constant independent of frequency. Also, in this model, a 
noise power spectra is assumed having amplitude o z , also Independent of 
spatial frequency. Linear equations resulting from this model do not 
require the existence of a well-defined Fourier transform and is as 
defined in Appendix A. There are advantages and disadvantages of this 
context-free approach. The major advantages are: 

• The theory permits a matrix multiplication approach with a fixed 
set of coefficients dependent only upon the antenna scan angle, 

• The method does not require any additional a priori knowledge, 

• The theory permits an adjustment of derived coefficients when the 
multiplication matrix is truncated. 

In fact, the theory demands that the sum of a truncated set of coeffi- 
cients should add up to one, and that the adjustment should be made addi- 
tively to each of the coefficients in the array. 
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There are also a number of disadvantages to a context-free process: 

• The theory is applicable to a single set of statistical assump- 
tions so that a theory which might work well in open ocean will 
not be optimal in a situation where there may exist water/land 
boundary. Hence, water/land interface artifacts will be present. 

• Optimization of the process to minimize land/water interface 
artifacts will sacrifice resolution in open ocean. 

• No context-free theory can improve spatial resolution without 
simultaneously incurring the disadvantages of a point source 
beaming energy into one of the side lobes of the resultant 
antenna point spread function. This is the Gibb's phenomena 
associated with the ringing that occurs near boundaries and is 
unavoidable in the region of any brightness temperature discon- 
tinuity. 

One of the major results of our previous analysis on context-free 
approaches is that avoidance of artifacts resulting from the antenna 
pattern correction process near land/water interfaces would require a 
point spread function which is not too different from the assumed model 
antenna point spread function. As such, a case is made for not 
performing any antenna pattern correction at all if one desire;, to avoid 
ringing artifacts. Applying a quadratic penalty function in the 
Backus-Gilbert cost function also leads to similar results. 

It therefore appears that for a context-free theory we either put up 
with ringing artifacts near land/water boundaries and take advantage of 
higher spatial resolution in mid-ocean, or, we avoid artifacts caused by 
the antenna correction process and put up with a less than optimal spa- 
tial resolution over relatively uniform areas. The best of both worlds 
can only be achieved with a context-sensitive antenna pattern correction 
process where the influences of land in the side lobes are compensated. 
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Context -Sensiti ve Theory 

A context-sensitive theory incorporates a priori knowledge of the 
real world into the antenna correction process. The coefficients to be 
utilized in an antenna pattern correction process could be adaptive based 
upon observed patterns in the data but more fundamentally, we propose a 
land/water geographic model be used to compute an a-priori antenna tem- 
perature world map. Fundamentally, a context-sensitive theory will 
pre-examine the array of antenna temperatures locally, and classify the 
scene to be one of several anticipated patterns based upon our knowledge 
of the real world. Once the scene is classified, a specific set of 
coefficients for antenna pattern corrections of any unexplained antenna 
temperature residuals will be applied which match the statistical charac- 
ter of the classified region. In such a theory, it may not be necessary 
to have a detailed knowledge of the locations of land/water interfaces 
since, to a large extent, these locations could be extracted from the 
observed behavior of antenna temperatures as they vary spatially. 
Nonetheless, it may be computationally more effective to utilize the 
geographic location data because no special effort would be needed for 
pattern recognition. 

In all of our context-sensitive approaches, a spatial model of 
brightness temperatures is required. The brightness temperature model 
may be obtained in alternative ways: 

• Land/water flags can be used from an earth surface map with an 
associated latitude-dependent brightness temperature model for 
both water and land. 

• An a priori brightness temperature model can be created from the 
LAMMR antenna temperature measurements which are processed and 
statistically incorporated in a world brightness map. 
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• Brightness temperatures for water and for land can be inferred 
from the antenna temperature measurements on the basis of land/ 
water flags and thereby utilized in a context-sensitive a priori 
parameterized model. 

• A bi -modal brightness temperature model can be assumed everywhere 
and utilized to avoid Gibb's artifacts that occur near land/water 
boundaries as well as at other geophysical fronts which may occur 
in nature. 

In each of the four aforementioned context-sensitive models, the 
approach to antenna pattern corrections would be to infer a priori 
antenna measurements from the assumed brightness temperature model. 

These modeled antenna temperatures are then subtracted from the observed 
antenna temperature measurements, yielding a variational antenna temper- 
ature measurement which is not explained by the assumed model. These 
unexplained antenna temperature fluctuations are then processed in a 
context-free spirit, utilizing a philosophy which is statistically 
appropriate . the prevailing conditions. One could expect that the 
observed fluctuations should have no a priori correlation and hence, a 
conservative minimum variance philosophy may be appropriate for enhancing 
the spatial degradations imposed by the modulation transfer function of 
the antenna. The inferred brightness temperature fluctuations can then 
be added to the brightness temperature a priori model to yield our best 
estimate of observed brightness temperatures. It is to be noted that 
many benefits can be obtained utilizing a context-sensitive approach to 
antenna pattern corrections. Some of the benefits are: 

• Discontinuities in the a priori brightness temperature model 
introduce high frequencies Into the interpretation of the antenna 
temperature measurements in a manner reminiscent of Shell's 
introduction of positive definite constraints in the interpreta- 
tion of radio telescope data, thereby avoiding the ringing arti- 
facts due to Gibbs' phenomena. 
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• Because an antenna pattern correction process Is applied on the 
fl actuations in antenna temperatures, we can relax our concerns a 
bit on the statistical assumptions that have been involved in 
converting these fluctuations to brightness temperature's fluctu- 
ations. 

• The theory could make use of a world brightness map, thereby 
justifying an approach which has many other side benefits espe- 
cially in the area of Quality Control and Assurance. 

• Such a theory directly attacks the problem of side lobe inter- 
ferences due to land masses, (side lobe compensation) and thereby 
permits the Isolation of the problem of spatial resolution 
enhancement which can now become a simplified process. We have 
shown it can be handled adequately utilizing a 5*5 or even a 3*3 
matrix multiplication approach. (It is also viable not do 
perform any spatial enhancement of the antenna temperature 
fluctuations at all, thereby significantly reducing real time 
computational burdens of the antenna pattern correction process.) 

Thus far, we have addressed our concerns on artifacts that are gen- 
erated due to the sharp spatial variation of brightness temperatures over 
the surface of the earth which the antenna does not adequately respond 
to. However, there are other sources of antenna temperature misinterpre- 
tations caused by biases which may be introduced into the system from the 
long tails of the point spread function response which although is low in 
amplitude, does however, cover large spatial extent. The integrated 
effects of the antenna temperatures caused by the tails of the distribu- 
tion are not corrected using a matrix approach to antenna pattern correc- 
tions. Experience with SMMR data has indicated such biases are related 
to the direction of coast lines relative to the flight path of the 
spacecraft. An a priori world brightness map could be utilized to 
compute biases caused by the tails of the antenna pattern, and these 
biases would then be Incorporated In the a priori estimate of antenna 
temperatures from the world brightness map, and Incorporated into the 
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antenna pattern correction process. Details of the four aforementioned 
approaches to context-sensitive processes are shown in Appendix A, 

Section A. 11. Analysis of these four approaches was conducted during the 
course of this work and the following paragraphs and comments are 
pertinent: 

Latitude Dependent Model 

The latitude dependent model establishes brightness temperatures for 
various catagorles of ocean, ice, snow and land which dependents both on 
latitude and seasonal (time). The use of such a model, of course, would 
only be approximate and local conditions would occur in which the actual 
brightness temperature, for example over land, might be significantly 
different than the global model. Under these circumstances the artifacts 
introduced by the land water discontinuity would not be total compensated 
for by our context sensitive model. The advantages, however, of the 
latitude dependent model is that it is parameter free and therefore not 
subject to errors that might occur in the datastream as might occur in 
any parameterized model. 

Parametric Land/Water Model 


In a parameterized model, the values of brightness temperatures are 
Inferred from the actual antenna measurements over some local array. The 
advantages of this model are that the parameterized values are extracted 
directly from the datastream and hence the discontinuity of brightness 
tempeatures across water, land boundaries are more accurately 
represented. The disadvantage of this scheme is that some degree of com- 
putation are required in tracking the brightness temperature parameters 
over the world surface types. Methods of establishing these parameter- 
ized brightness temperatures are shown in Appendix A, Section A. 11. 

By -Modal Model 

In the by -modal mode! no a priori knowledge of the world structure 
is required. The algorithum continuously assumes that each region is 
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by-modal consisting of two species. A threshold temperature is estab- 
lished to distinguish both species. Brightness temperature estimates are 
established for both species and an antenna temperature is therein 
computed. The disadvantage to this scheme is that it is computationally 
complex and not really applicable over open ocean or open land. It has 
however, the advantage of representing weather front discontinuities 
which the other techniques cannot accommodate. (Large weather fronts, 
however, are not really modelable as discontinuities in brightness 
temperatures so that this later advantage is not truly applicable). 

World Brightness Map 

In the world brightness map approach the a priori model of 
brightness temperature is obtained from an analysis of data itself and 
stored in a statistical data set. The antenna temperature expected from 
a given region is extracted and subtracted from the observed field. This 
approach suffers from the disadvantage that statistical fluctuations in 
brightness temperature over land would still yield discontinuity 
artifacts and also from the disadvantage that the process requires a 
continuous updating from a world brightness map and associated 
statistics. It has, however, the advantage that such a world brightness 
map would have other benefits to the LAMMR ground processing systems. 

A tradeoff analysis was conducted for each of these four systems 
from the expected error point-cf-view. It was determined that the para- 
metric land/water model was most appropriate for the LAMMR processing 
system. 

3.2 Process Overview 


Figure 3.0 shows a block diagram of the context sensitive antenna 
pattern correction and side lobe compensation processes. Details of the 
APC and SLC process is shown in Figure 3.1. Here the a priori brightness 
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temperature map is shown being convolved with the antenna point spread 
function to generate the a priori antenna temperature (step 2). The 
measured antenna temperatures of the scene Is modified by subtracting the 
expected antenna temperature from the model generating an antenna 
temperature variation. These three steps constitute the side lobe 
compensation process. In that, a priori knowledge of the Earth's geometry 
is Incorporated so as to remove the side lobe influences of land. In 
step 4, brightness temperature variations are obtained through an antenna 
pattern correction process of the matrix variety and these are added in 
step 5 to the a priori brightness model generating an estimated 
brightness temperature. These brightness temperatures are self 
consistent, in that, a convolution with the antenna profile would 
generate the original antenna temperature as observed. An alternative 
approach to context-sensitive antenna pattern corrections is shown in 
Figure 3.2. It uses different APC processes in different circumstances. 
Local regions are examined to determine their homogeneities by accessing 
information from the a priori Spatial Model. Spatial frequency 
enhancement processing is applied to homogeneous regions (ocean) 
side-lobe suppression processing is applied to heterogeneous regions 
where there are large signal fluctuations (land-ocean boundaries). 

In studying the details of this approach, it was originally 
conceived that different APC coefficients would be utilized in a 
relatively few sets of circumstances, such as over ocean, over land a, id 
at land-ocean boundaries. However, it was soon realized that proper 
consideration at the land ocean boundaries would require a further 
disection according to boundary directions. Although the approach is 
viable, the large number of antenna pattern correction coefficients that 
would be required to handle the different spatial circumstances would be 
overwhelmingly large so that this alternative approach is not 
recoranended. 
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Figure 3.2 APC Context-Sensitive Process 
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3.3 Generation of Surface Type Map 


The generation of surf, a type map (see Figure 3.3) involves the 
segmentation of various geophysical features and land masses from the 
ocean. It car be done by thresholding and utilizing known geographical 
data. An alternative approach to the generation of a surface feature map 
is to utilize surface feature boundary tapes and convert them to a grid 
of latitute/longitude surface types. A code exist at Goddard Space 
Flight Center for performing such a function. However, as previously 
mentioned, such a conversion is non-trivia! and to date the conversion 
process leads to rather wide dispersion near the boundaries. 

Nonetheless, it is believed that modifications of this code is possible 
with relatively little effort and presents perhaps a better alternative 
to the thresholding procedure described above. 

3.4 Our Planning Effort 


Figure 3.4 shows a summary of our antenna pattern correction plan. 
It consists of five basic areas. 

• Context-free antenna pattern correction 

• Context-sensitive antenna pattern correction 

• Synthetic scene simulation 

• Performance evaluations 

• Documentation 

A brief description of this plan follows: 

APC Processes 


- correct side-lobe effects (SLC) 

- enhance spatial frequencies (APC) 
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(A) Context-free APC 

- Context-free APC algorithms 

- Performance evaluation 

- Prototype software 

- Sets of correction coefficients 

(B) Context-sensitive APC 

- Context-sensitive APC algorithms 

- Sets of correction coefficients 

- Performance evaluation 

- Prototype software 

(C) Synthetic Scene Simulation 

- Various simulated scenes 

- simulator software 

(D) Performance Evaluations 

- Context-free APC apply to simulated scenes 

- Context-sensitive APC apply to simulated scenes 

- Comparison/evaluation of APC responses on simulated scenes 

- Context-free APC apply to SMMR data 

- Context-sensitive APC apply to SMMR data 

- Comparison/evaluation of APC responses of SMMR data 

- Trade-off evaluation of context-free vs context-sensitive ao^roaches 

(E) Docunentations 


- Final report and presentation 

- Monthly progress reports 

- Documentation of all software nodules 

Figure 3.4 Summary of APC Plan 
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(i ) Context-free APC 


A linear-filter APC algorithm based on least-square 
constraints has been developed and implemented. It enhances 
the spatial responses of the antenna pattern and at the same 
time ehances the Gibb's oscillations near discontinuities 
(example: land-ocean boundaries). This APC philosophy 

employs no geographical information and hence suffers from 
the fact that the antenna side lobes can be suppressed only 
at the expense of spatial resolution. 

( i i ) Context-sensitive APC 


It is feasible to suppress the antenna side lobes as 
well as enhance spatial frequencies by incorporating some 
a priori information into the correction process. One way 
to achieve this is to utilize information about where land 
masses are situated geographically and approach the 
correction problem in a pattern recognition point of view. 

This involves the using of a-priori models of antenna 
temperatures and comparing these with observation. 

Context-sensitive APC Development 

(i) Perform analysis for context-sensitive APC processes 

(ii) Define scenarios on how the context-sensitive processes are to be 
performed 

(iii) Context-sensitive processes error evaluation and selection 

(iv) Generate context-sensitive APC coefficients 

(v) Software development of the APC process 

- design 

- coding 

- testing 

- documentation 
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Antenna Temperacure Simulator 

(i) Generate a synthetic scene with land/ocean boundaries and 
various geophysical parameters. 

(ii) Generate Tg's that correspond to the given geophysical 
parameters. 

(iii) convert the given scene of Tg's into T^'s with the 
effects of the side-lobes of the antenna. 

(iv) Software development of the simulator 

- design 

- coding 

- testing 

- documentation 
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Testing and Evaluation 

The software will take as input both a simulated scene and the SMMR 
data. Both the context-free approach and the context-sensitive approach 
will perform the APC to retrieve the Tg's. A trade-off analysis will 
be performed to evaluate the two APC approaches based on the responses of 
the simulated scene and of the SMMR data. 



Work Performed 


Most of the effort progressed on schedule as planned, except for the 
use of SMMR data in replacing our simulated scene. At this juncture, we 
had difficulties obtaining a useable surface type map which delineated 
land water boundaries sharply and also discovered that the SMMR data, 
because of the binning process, no longer revealed the point spread 
function characteristics of the antenna; it would not constitute a valid 
test for evaluation. Hence, more effort was spent on the simulated scene 
approach. 
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As a result of the evaluation, we have shown that the context- 
sensitive antenna pattern correction process as originally conceived 
indeed avoids many of the aliasing artifacts that are generated in the 
context-free antenna pattern correction process. Hence, it is highly 
recommended that LAMMR processing include the use of a brightness 
temperature world map for side-lobe compensation and subsequent antenna 
pattern correction for enhancing spatial frequencies. 

3.5 Analytic and Simulated Results 

The context-sensitive side lobe compensation and antenna pattern 
correction process was simulated on a synthetic image shown in Figure 3.5 
called the rectangular island. In this figure, the brightness 
temperature over water was designated by a zero while the brightness 
temperature over land was designated by a numeric 4. These values were 
chosen so that the island could be represented as a single digit output. 
Mapping these values into realistic brightness temperatures is easily 
accomplished by applying a gain and an offset to the rectangular island 
file within the computer. The simulation plan calls for comparing the 
context-free antenna pattern correction process with a context-sensitive 
one. In achieving this end, it was decided that the k-space Fourier 
domain would be the most suitable space within which to do the required 
computations. The fast Fourier transformed algorithm was utilized, and 
Figure 3.6 gives a representation of the rectangular island in the 
k-space Fourier domain. Fourier transforms were taken with and without 
additive noise to the simulated antenna temperature field. This field 
was obtained by multiplying the ideal antenna transfer function by the 
Fourier transform of the rectangular island image and performing the 
inverse Fourier transform. Figure 3.7 shows a trace of the antenna 
temperature across the boundary of the island. The antenna was assumed 
to be a circular aperture with no tapering thereby leading to a first 
order Bessel function type point spread function as described in the 
previous report. Note how the antenna function blurs the discontinuity 
between water and land and that the residual side lobes cause some minor 
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oscillation as the antenna temperature approaches the steady state values 
of land and water. A 5*5 context-free antenna correction process was 
applied to an array of extracted antenna temperatures obtained selecting 
every fourth value for every fourth line. This sampling strategy 
corresponded to a minor amount of undersampling as in the design of the 
LAMMR Instrument. Figure 3.8 shows the trace of the context-free antenna 
pattern correction process applied to the antenna temperatures for 
estimating brightness temperatures. Prior to the antenna correction 
process a realistic amount of noise was added to the antenna 
temperature. From the results we can see traces of the Gibb's phenomenon 
which in this example is a bit subdued due to the slight undersampling 
and therefore, inability to totally compensate frequencies up to the 
antenna cutoff frequency. The context-sensitive algorithm which has been 
previously described was applied to the same noisy array of antenna 
temperatures for the rectangular Island. A trace across the island water 
boundary interface is shown in Figure 3.9. As expected, the simulation 
shows a reconstruction of the original brightness temperature scene with 
minor perturbations most probably induced by the additive noise. Many 
experiments were conducted, the behavior of which could be totally 
explained on the basis of previously developed theories. In general the 
results are not too startling and in some respects, many of the potential 
advantages of a context-sensitive implementation for brightness 
temperature estimates have not really been tested. In addition to 5*5 
matrices, 7*7 and 3*3 arrays sizes have been tried with minor differences 
of the ones already presented. Figure 3.10 shows the frequency response 
of the antenna pattern correction arrays associated with different matrix 
sizes. In the higher frequency region, one can see that 5*5 array 
behaves as well as a 7*7 and that the 3*3 antenna pattern correction 
compensates for high spatial frequencies in a manner just slightly below 
the optimum nominal. 

In conclusion, the context-sensitive antenna pattern correction and 
side lobe compensation processes herein described generate results which 
are generally superior over the context-free approaches, while simultane- 
ously reducing the amount of real time computation albiet the need to 
generate a brightness temperature world data base. 


49 




LAND 


Business ami Technological Systems , Inc 



50 


Figure 3.7 Trace of the Antenna Temperature 
Across the Rectangular Island Boundary 
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Figure 3.8 Trace of Context-Free Antenna Correction Process APC 
For Estimating Brightness Temperatures (5x5) with Additive 

Noise 



Figure 3.9 Trace of ^text-Sen^itive Antenna Correction Process For 
Estimating Brightness Temperatures (5x5 APC on <5T A ) with Additive Noise 


magnitude 


Bisixtss i \i) Tu hsoum.hal Systems, Im: 



Figure 3.10 Frequency Response of Antenna 
Patterns Corrected by Different Matrix Sizes 
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4.0 SOFTWARE DESCRIPTION 


In this chapter, the implementation of APC software and the demon- 
stration of context dependent APC philosophies are described. 

In prior work, context-free philosophies were implemented and 
evaluated. A model of two-dimensional power pattern was generated upon 
which the APC process was performed. A detailed geometric transformation 
was developed and coded. This transformation relates the antenna pattern 
as might be provided through antenna testing to the instrument coordinate 
system in which digital data is sampled. A system of software was 
developed to generate antenna pattern correction coefficients in 
conformance with the established context-free analysis. The system 
includes routines to calculate cross-correlation, to shift and rotate 
antenna patterns, and to solve linear equations. Software to generate 
effective point spread function in instrument coordinates was also 
developed. Another set of routines was implemented to perform Fourier 
analysis of the correction processes. With some enhancements, this 
system provides the basic tool to demonstrate and evaluate the proposed 
context-dependent antenna pattern correction and side lobe compensation 
philosophies. 

The design and implementation of the APC software functions are 
based on top-down and modular design. This structural desiqn will be 
discussed in detail in Section 4.2. The code used in this study is given 
in Appendix B. 

4.1 Functional Organization and Overview 

Figure 4.1 shows a system overview of the performed functions. 

The programs were written in FORTRAN and implemented on the IBM 360/91 at 
NASA/GSFC. The operational software consists of the following functions: 
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Figure 4.1 Software System Overview 
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(1) Simulation of an antenna power pattern on a 64*64 plane 
(k-space) parallel to the antenna aperture. 

(2) Transformation of coordinates from the k-space to the instru- 
ment coordinates (x-space). 

(3) Generation of necessary parameters required by the optimization 
processes. These parameters are used in the data interpretive 
equations and the matrix truncation approach. 

(4) Determination of correction coefficients by solving the set of 
linear equations. 

(5) Correction of antenna patterns to generate various point spread 
funcvlonsx 

(6) Perform two-dimensional discrete Fourier transform of various 
point spread functions. 

(7) Generate synthetic image 

(8) Two-Dimensional Fourier Transform on Synthetic Image 

(9) Generate antenna temperature field 

(10) Add random white noise 

(11) Perform matrix HPC operation 

(12) Difference simulated T^ and <T^> model 

(13) Perform matrix APC operator 

(14) Compute brightness temperatures 

(15) Generate synthetic scene 


In the following subsections the above-mentioned functions are 
described in detail and flow diagrams are provided. 


4.1.1 Simulation of Antenna Pattern 


Using the concepts of linearity, the impulse response of an incoher- 
ent system due to a circular aperture with a 0-1 function having 1 in the 
area of the aperture and 0 elsewhere is the so-called Airy disk. In 
classical optics, an Airy disk is the diffraction pattern intensity for a 
uniformly illuminated circular aperture. The analytic solution to this 
pattern is the square of a first order Bessel function divided by its 
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argument, [Jj(k)/k] . By assuming a 0-1 real function circular antenna 

aperture, an antenna power pattern formed on a plane parallel to the 
antenna aperture will have an Airy disk characteristic. 

Subroutine LENS was implemented to generate a 64*64 Airy disk which 

2 

corresponds to [J^(k)/k] . The diameter of the disk was chosen to be 32 
with 64 equally-spaced intervals along its diameter. The disk contains 
the main lobe and more than four side lobes as shown in Figure 4.2(d). 

The power at each location is computed by the polynomial approximations 
as shown below: 


[Jl(k)l - I k 2 It 4 

-y— ' * 7 - 0.56249 (y) + 0.21093 (y) 

- 0.03954 (y ) 6 + 0.00443 (y ) 8 

,jt 10 , k 12 

- 0.00031 (y) + 0.00001 (y) 


for 


-3 < k < 3 


wnere 


k 



+ k y 


2 


with 


For 


where 


k k = -15.5, -15.0, -14.5, •••0.0,»««, 15.5, 16.0 
x y 


3 < k < ® 


[Jj(k)3 ^12 

Jy— * k f Cos 0 


3 3 2 

f =0.79788 + 0.000001 (y) + 0.1659 (y) 

+ 0.00017 (|) 3 - 0.00249 (y)* 
+ 0.00113 (y ) 5 - 0.00020 (y ) 6 
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Figure 4.2 Airy disk 


a. 

b. 

c. 

d. 


Circular antenna operture 

J. (k) 2 

Antenna power pattern, [-^ — ] 

Perspective uf the antenna power pattern 
One-dimensional plot of the disk pattern along its diameter 

J,(x) 2 

[— f 

L x J 
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and 


3 3 2 

9 « k -2.35619 + 0.12499 {^) + 0.00005 (|) 

3 3 34 

- 0.00637 {-£) - 0.00074 (|) 

+ 0.00079 (|) 5 - 0.00029 fy* 

The generated antenna power pattern was normalized to have a total power 
of one and the resulting Airy disk was stored in a disk file named 
AIRY. DATA using 1/0 unit number 12 for ready access. The data flow is 
shown in Figure 4.3. 


Radius of 
disk 



Figure 4.3 Data Flow of Subroutine LENS 

TESTLS: Test program for subroutine LENS 

LENS: Routine to generate the Airy disk 
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4.1.2 Coordinate Transformation 


Our antenna power pattern is presented in a plane (k-space) 
orthogonal to the pointing axis of the antenna. We need to transform 
this coordinate system to instrument coordinates in which the grid is 
represented by a differential change along the scan and across the scan 
line. 


We begin with the scan pixel i and line number j in instrument 
coordinates and compute the antenna power and its location in the k-space 
by the following transformation procedures: 

(1) The antenna power pattern (the Airy disk) is retrieved. 

It consists of a two-dimensional array 

G(k x , k y ) 

where (k , k ) i; the integer grid point and G is the 

X J 

antenna power at (k , k ) . 

x y 

(2) Given Sg , the unit vector of the boresight direction of the 
antenna, the position of the antenna with respect to the 
earth, and the instrument coordinates ( i , j ) , the unit 
vectors S from the antenna to (i,j) can be computed. 

(3) Then, a pointing vector £ projected in the k-space ortho- 
gonal to Sg is expressed as 

R • [(S-S 0 ) - (S-Sq) • $ 0 S 0 ] where 

D is the antenna operture diameter and X the wavelength. 
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(4) From k, we can compute the x and y components of k which 

permits us to Interpolate the antenna power pattern G(k , k ) 

x y 

where (k ,k ) Is closed to (x,y). 

X Jr 

(5) The relationship between the power pattern in the k-space and 
the power pattern in the instrument coordinates Is given by 

A ^ 

A = G k 2 Cosv 0 v o r Q AtA0 Cos(e-0 Q ) 

S 

where A is the power pattern in instrument coordinates, 

«q is the ground velocity of the satellite subtrack which depends upon 
earth rotation. It is latitude dependent and assumed constant in the APC 
computations (Vg = 7 km/sec.) 

rg is a nominally constant radius of the ground trace scanned circle. It 
depends upon the distance from the antenna to the earth's center, the 
earth's radius and the instrument cone angle ^g. 

At is the time between adjacent line C At = 1 sec) 

A0 is the separation in radians between adjacent pixel samples 

0g is the angle at which the sweep is orthogonal to the ground velocity 
vector. In the simplified model, we assume there is no earth rotation 
and hence 0g a 0. 

s is the Instantaneous pointing unit vector from the antenna to the pixel 
sample. 

p Is the instantaneous pointing unit vector from the earth's center to 
the pixel sample. 
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£ Is the pointing vector projected onto the k-space 

<J*q Is the Instrument pointing angle from Nadir ( <|»q = 43.6°), and 

G is the power pattern in k-space. 

Subroutine TXFORM was Implemented to transform an instrument 
coordinates (1,j) to (x,y). Subroutine INT3P was implemented to locate 
three integer grid points in the k-space which are closest to (x,y). 
Subroutine INTLN performs the bi-1 inear Interpolation to obtain an 
interpolated power magnitude at (x,y). The main program ANTENA was 
Implemented to link all three subroutines together to generate a file of 
64 x 64 elements containing an antenna oower pattern in instrument 
coordinates. The grid spacing Is 7/4 kilometers. The generated pattern 
named ANT. DATA was normalized to have a total power of one and was stored 
in a disk file using I/O unit number 11. Figure 4.4 shows the data flow. 



Figure 4.4 

Flow diagram of the coordinate transformation and the generation of the 
antenna power pattern, ANT. DATA, in Instrument coordinates. 
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(1,j) Is the Instrument coordinates 
(x,y) Is the k-space coordinates 

* * 

constants * k 2 Cos 4 » 0 v 0 r 0 AtAQ Cos(o-9 Q ) — 

G * power pattern In k-space 
A * power pattern In Instrument coordinates 
INT3P and INTLN are routines to Interpolate A 


4,1.3 Data Interpretive Equations and Matrix Truncation 

A modeled approach Is u*_d to estimate the brightness temperature. 

2 

It Involves the estimation of the optimum set of coefficients M that 
minimizes the least squares error in the temperature fluctuations. 

Taking a functional variation of the variance with respect to the 
estimation coefficients leads to a set of normal equations. By assuming 
the statistical independence of the brightness temperature fluctuations 
and a delta function burring effects of the antenna, we simplify the set 
of normal equations to an expression that only Involves the cross- 
correlation function and shifting of the antenna pattern. We further 
select a finite set of grid points and truncate the grid, the simplified 
normal equations become a simple matrix equation: 

= == 2 = -1 
M = r(£ + a I) 


where M is the estimation coefficients, 

3 

Z represents the cross-correlation function cr the antenna pattern 
2 

iMs the shifted antenna pattern, and 
a Is the signal-independent noise. 


An element Zj j of the matrix z is approximated by 

64 64 


c u* s ‘ Ji J/'V.’ Vn> W 


An Ant 
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where i * (i , i ) and j = (j , j ) are locations of sample pixels, (n,m) 
x y x y 

is an Integer grid point of the antenna pattern, An, Am are spaclngs of 

1 2 

the grids (An * Am = -^1, S is the signal of the brightness temperature 
fluctuation, and A is the antenna power pattern. In estimating the 3x3 
coefficient matrix, 2 Is a 9x9 matrix with elements arranged in a fashion 
as shown in Figure 4.5. 

3 

A 3x3 coefficient matrix M 


A 
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3 


The 9x9 cross-correlation matrix 2 
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Figure 4.5 

3 

The 2 matrix of a 3x3 matrix equation 
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The elements of the r matrix is nothing more than the shifted 

o 

antenna power pattern multiplies by S' given by 

r i - s 2 A <V V 

In the 3x3 coefficient matrix problem, r is a 9x1 column vector and can 
be written as ? = ,r i^‘ 

The matrix equation then becomes 

m = r *s 2 (rs 2 + o 2 I )’ 1 



s p 3 = 2 z <j2 

where E'S * £, P'S - r and is the square of the noise to signal 

S « = 

ratio. This equation is in the form of AX * B linear equations where 

given the A and 1 matrices, 1 can be solved for by Gaussian elimination 
with partial pivoting. 

Subroutine CONVOL was implemented to compute the cross-correlation 
of two antenna patterns centered at i and j respectively. Both i and j 
are neighborhood sample pixels of the center of a 64x64 plane. The 
computed cross-correlation is used to generate the S' of various sizes. 

x 

The 2' matrix can be a 9x9, 25x25, or 49x49 matrix determined by the size 
of the matrix of correction coefficients. Figure 4.6 shows three 
different window sizes of M with respect to the antenna pattern. 

x 

Subroutine GAMMA was implemented to compute the matrix r*. The 

x = 

generated P' is a 9x1, 25x1, or 49x1 row vectors. Having obtained £' and 

x 

r', subroutine COEF solve the linear equations by using the IMSL 

a 

subroutine LEQT1F. The resulting M contains correction coefficients and 
they are normalized to have a sun of one. 
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@ center of a 64x64 antenna pattern 
O sample pixel s 

grids of antenna pattern 


Figure 4.6 

The 3x3, 5x5, and 7x7 correction windows superimpose on a 64x64 antenna 
pattern. 
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Figure 4.7 presents three sets of coefficients generated by the above 
described procedures with 10.65 GHZ and a signal to noise ratio of 10. 
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Figure 4.7 

Correction Coefficient Matrices M, of sizes 3x3, 5x5, and 7x7 at 
u> » 10.65 GHZ and (o/S) 2 * 0.01 
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4.1.4 Point Spread Function of the Corrected Antenna Pattern 

The optimum set of correction coefficients R can be used to estimate 
the brightness temperature measurements Tg by the linear combination of 
R and the antenna temperature measurements T ft . It has been recognized 
that the estimated Tg is associated with a specific point spread function 
given by 


PSF(x) = I Mr A(x + 1) 

i 1 

where R can be thought of as a local template that operates on each and 
every pixel x of the antenna pattern A. The vector 1 points from the 
center of the template to pixels inside the template and M-. represents 
the coefficient at i. The value of the PSF at x is then computed by the 
weighted sum of A at various locations (x + 1). 

Subroutine ANTPSF was implemented to comput the point spread 
function of the antenna pattern ANT. DATA. For the 3x3 coefficient 
matrix, 37,000 multiplications and additions are needed to generate the 
64x64 PSF . The resultant point spread function is stored on a disk file 
named PSF. DATA using I/O unit number 13. Figure 4.8 shows the data flow 
and all the software modules required to generate the coefficients and 
point spread function. 

A system of software was implemented to perform Fourier analysis of 
the antenna pattern correction processes. It includes routines to 
perform two-dimensional discrete Fourier transform, and to compute 
decibels, power spectrum and various output functions. 

4.? Software Structure 

Top-down design and modular programming are used in the design and 
implementation of the APC demonstration software. All the routines are 
documented and listed in Appendix B. 
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SIZE * 3, 5, or 7 

< J 

— - noise to signal ratio 

(i,j) = pixel locations 

s 

Z = cross-correlation matrix 

s 

r = shifted antenna matrix 

2 

M = correction coefficients 

ANT =* antenna power pattern 

PSF = point spread function 

Figure 4.8 


Data Flow and Software Modules for the Generation 
of Coefficients and Point Spread Function 
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The idea of top-down design can be viewed as breaking down the 
problem into simpler problems and checking the relationships between them 
before proceeding to a more detailed specification. In programming, it 
involves decomposing the overall complicated program into precisely 
specified subprograms, so that each subprogram is solved correctly, and 
fitting these solutions together in a specified way. This programming 
method allows us to generate clean algorithms with far fewer mistakes 
than by the conventional line-by-line method. 

Modular programming has several objectives. (1) A program module 
must be working correctly by itself regardless of the context of its use 
in building larger systems of software. (2) To realize large software 
systems, program modules written under different authorities must be 
conveniently fitted together without changes in any of the component 
modules. A module is joined to other modules by communication links. 

Each module receives data from its input links, transforms it, and 
outputs it to other modules over its output links. Function subprograms 
and subroutine subprograms can serve as program modules. The use of 
FORTRAN is a reasonable choice in the implementation of the APC software 
to test out module design. 

Top-down modular analysis can be described by a graphical model in 
the form of a tree of the functions to be performed. Attached to each 
node of the tree are programs to perform specialized operations and 
structure tests. Each node can branch to one of the other nodes in the 
same level. All nodes are also modular, meaning they are completely 
independent of the other nodes. Programs that are designed on the basis 
of this method are very simple to debug or to modify. Each module can be 
tested separately and when it works, can be added to the system. 

The tree diagram of the software system for the APC processes is 
shown in Figure 4.9. This graphical structure specifies the types and 
instances when primitive operators in terms of subroutines are used 
during the APC processes. The program modules below the first level are 
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Figure 4.9 Software System and its Structure for the APC Demonstration 
Procedures (Descriptions of Software Modules are Listed In 
Section 4.2) 
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primitive operators that perform function described in previous 
sections. Attached to each module of the software tree structure are a 
list of program variables and common storage areas. All comnuni cation 
between modules is through the program variables connecting them at a 
higher or lower level in the tree. The modules at all levels are 
basically independent of each other. Hence, changes in any of the 
modules may only require modifications in themselves and a minimum number 
of changes in the nodes directly connected to them. The parent node is a 
supervising program. It performs the functions of linking the primitive 
modules together, branching to appropriate ones, testing their program 
variables, evaluating their performances, and outputting results. 

The following summarizes the various primitive modules implemented. 
Each is associated with a different primitive operation, and each is 
stored in a disk file. 


PRIMITIVE MODULE DESCRIPTION 


LENS To generate a 64x64 Airy disk, which corresponds 

to the square of the first-order Bessel function 
divided by its argument. 

INPUTS: RADIUS (radius of the disk) 

OUTPUTS: AIRY (airy disk on disk file) 

TXFORM To transform a discrete location in the instrument 
coordinates to an antenna coordinate system 
orthogonal to the antenna pointing axis. 


INPUTS : 


OUTPUTS: 


(FJX, FJY) (Antenna Pointing Direction) 
(FIX, FIY) (pixel location on ground) 

GHZ (radiometer frequency) 

NOSPS (number of samples/scan line) 
(QX, QY) (transformed coordinates) 

CONST (transformed antenna constant) 
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INT3P 


INTLN 


CONVOL 


SIGMA 


GAMMA 


To locate 3 nearest neighborhood grid points 

INPUTS : (QX, QY) (input point) 

OUTPUTS : (AX, AY) 

(BX, BY) (3 nearest grid points) 

(CX, CY) 


To perform bi-1 inear interpolation 


INPUTS : (QX, QY) 
(AX, AY) 
(BX, BY) 
(CX, CY) 
AIRY 

OUTPUTS : VALUE 


(input point) 

(3 nearest grid points) 

(Common block contains Airy disk) 
(interpolated value at (QX, QY)) 


To compute the cross-correlation of two antenna power 
patterns centered at i and j respectively. 

INPUTS : (1X1, IY1) (1st point) 

(1X2, IY2) (2nd point) 

OUTPUTS : El 1 12 (the cross-correlation) 


To compute the E matrix 

INPUTS : ISIZE (window size of M) 

ANT (common block contains power 

pattern) 

OUTPUT : EMATIX (matrix E) 

To compute the T matrix, a windowed antenna power 
pattern 
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INPUTS : ISIZE 
ANT 

OUTPUTS : GMATIX 

COEF To solve the set of 

INPUTS : ISIZE 
SNTSR 

EMATIX 

GMATIX 

OUTPUTS : GMATIX 

ANTPSF To compute the point 

specific coefficient 

INPUTS : ISIZE 
GMATIX 
ANT 

OUTPUTS : PSF 


(window size of M) 

(common block contains power 
pattern) 

x 

(matrix r) 

inear equations 

(window size of M) 

(square of noise to signal ratio) 

(common block contain E) 

(common block contain T) 

(common block contain M) 

spread function associated with a 

x 

matrix M 

(window size of M) 

(common block contains M) 

(common block contain power 
pattern) 

(common block contain corrected 
pattern) 


DFT2D 


To compute the two-dimensional discrete Fourier 
transform of an MXN image. 


INPUTS : (M,N) 
FXY 

OUTPUTS : FUVRL 
FUVIM 


(size of the input image) 
(input image) 

(transformed image, real) 
(transformed image, imaginary) 
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DFT2XY 


FSCENE 


To compute the two-dimensional discrete Fourier 
transform along the ? axes of an 64x64 image, to 
compute their power and dB, and to output results.. 


INPUTS : FXY 
OUTPUTS : FUREAL 
FUIMAG 

FVREAL 

FUIMAG 

FUMAGN 

FVMAGN 

FUDB 

FVDB 


(input image) 

(2-0 OFT real along V*0 axis) 
(2-D OFT imaginary along V=0 
axis) 

(2-0 DFT real along U=0 axis) 
(2-D DFT imaginary along U=0 
axis) 

(2-0 power spectrum along V=0 
axis) 

(2-D power spectrum along 11=0 
axis) 

(dB al ong V=0 axis) 

(dB along U=0 axis) 


To compute the fourier transform of an image 
using .,ubroute DFT2D. The resulting real and 
imaginary parts are stored at FSCENER. DATA and 
FSCENEI.DATA with 15 and 16 as LUN's 

INPUTS : IDUM (64*64 input image) 

OUTPUTS : Are Real and Imaginary 

transforms stored in LUN 15 
and 16. 
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MAIN PROGRAM DESCRIPTION 

TESTLN Test program for subroutine LENS and to create the 

Airy disk. 

ANTENA Main program to generate a simulated antenna power 

pattern on an instrument coordinate system. 

INPUTS : AIRY (common block contains Airy disk) 

GHZ (primary radiometer frequency) 

OUTPUTS : ANT (common block contains antenna 

power pattern) 


TESTSA Main program to generate the APC correction matrix R, 

to compute the corrected antenna point spread 
function, and to perform Fourier analysis. 


INPUTS 

: ANT 

(common block contains antenna 
power pattern) 


I SIZE 

(window size of correction matrix 
R) 


SNTSR 

(square of noise to signal ratio) 

OUTPUTS 

: PSF 

(common block contains point 
spread function) 
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PRECEDING PAGE BLANK NOT FILMED 


5.0 CONCLUSIONS AND RECOMMENDATIONS 

The analyses and demonstrations in this report illustrate a method 
for processing LAMMR data so as to extract a context sensitive 
brightness temperature for side lobe compensation and at the same time, 
to modify the antenna pattern for frequency enhancements. The objectives 
in this processing are to suppress the influences of the side lobes and 
to enhance the spatial frequency responses. The influences from side 
lobes become a severe problem when the antenna encounters land -ocean 
boundaries, since the contrast between land and ocean is very high. 

There is, however, a serious problem in connection with the antenna 
patterns: there is a loss in spatial responses as a result of the finite 

size of the receiving antenna. In this context, the ideal antenna 
pattern correction processing should correct the side-lobe influences and 
enhance the spatial frequency responses. Both of these desires cannot be 
simul taneously achieved with linear processing; context sensitive 
processing is the answer. The antenna side lobes can be suppressed by a 
linear filter only at the expense of spatial resolution. On the other 
hand, the enhancement of the pattern also enhances Gibb's oscillations 
near discontinuities. There is no way to compensate the antenna 
frequency responses and avoid radiance from nearby land masses beaming 
through the antenna side lobes. However, this study shows it is feasible 
to suppress the Gibb's oscillations as well as enhance frequences by 
incorporating an a-priori model in the correction process. This can be 
done by using information about where land masses are situated 
geographical ly. In so doing, a process has been established to remove 
side-lobe influences near land-ocean interface and to simultaneously 
enhance spatial frequency responses everywhere. The results of this 
study show how to use the a-priori knowledge of location of land-ocean 
boundaries for antenna pattern correction. 

It has been shown that no single context-free antenna pattern 
correction process will adequately compensate for the antenna side -lobe 
influences, while simultaneously enhance the spatial frequency responses 
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modulated by the physical limitations of the antenna. It is for this 
reason that a context-sensitive antenna pattern correction process 
was developed in details and was implemented to prove its feasibility. 
There are several context-sensitive techniques which were applied in this 
application study. 

• A statistically-derived world brightness map that contains a priori 
knowledge of expected brightness temperature over the surface of the 
earth was simulated to generate expected antenna temperatures and 
the temperature biases associated with the geographical locations of 
land and ocean. These biases were incorporated into the context- 
sensitive APC synthetic model to compensate for the side-lobe 
effects. 

• The use of a two-model procedure in conjunction with the world 
brightness map was simulated for the applications of different sets 
of coefficients in two different situation. In the case of near 
land-ocean boundaries, a special set of correction coefficients was 
applied so as to eliminate some of the artifacts associated with 
Gibb's oscillations. Another general set of correction coefficients 
was then be used everywhere over uniform ocean to enhance spatial 
frequencies. The set near the boundaries, however, blured that 
interface and did not really meet the objectives. 

• Another approach was to hypothesize the existence of land-ocean 
boundaries everywhere and apply a set of coefficients that 
suppresses side lobes to the entire scene, thereby affecting the APC 
performance over the ocean. Analysis proved this technique to be 
too complex for Real time ground processing. 

Errors measurements of these different context-sensitive approaches 
were assessed and a trade-off analysis was conducted to make the final 
selection. The trade-off included performance evaluations. Computa- 
tional requirements such as resolution versus side-lobe suppression and 
noise amplications were also studied and used in the trade-off 
eval uation. 
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The Recommended Algorithm 

In general, a model of expected brightness temperatures , <Tg>, must 
be constructed based on the geographical locations of land-ocean 
boundaries and the statistical data of ocean phenomena. Given this world 
brightness model, the expected antenna temperatures, <T A >, can be 
generated. The fluctuation of antenna temperature, 61^, are then 
computed as the difference between expected and measured antenna 
temperatures. APC process matched to the statistical data base are 
then be applied to the fluctuation 6T^ in order to generate the 
fluctuation brightness temperatures 5Tg. The final data interpreted 
brightness temperature is computed as: 

t b = <t b > + 6T b* 

It is recommended that this context-sensitive approaches be 
implemented for LAMMR ground processing. A complete evaluation of 
various context-sensitive approaches and context-free approaches showed 
it to be the most practicable. 
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Appendix A 
Detailed Analysis 

A.l Basic Mathematical Formulation of Antenna Temperatures 

In this section, we deal specifically with the problem of performing 
transformations on the antenna point spread functions that are viable in 
order to maximize the information contained in the data without being 
severely penalized by the influences of noise. 

There are several elements of concern in the 'ormance of the 
antenna pattern correction: 

• The desire to limit the convolution array size needed to 
perform the APC. 

• The formulation of the problem in the continuum domain so 
as to avoide complications that may result from aliasing 
inf! uertces. 

• The desire to include alternative measures of performance 
as suggested by A. Stogryn. 

In general , the problem that must be solved involves the inversion 
of a set of horizontal and vertical polarization antenna temperatures 
T^ (2x1 vector) that are sampled uniformly over some grid to an estimate 
of the brightness temperature components of polarization any* ere in the 
spatial fielo. In usual applications, the brightness temperature esti- 
mates are evaluated at the same set of grid points upon which the 
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measurements of antenna temperature are made. However, in general , the 
estimation problem also provides for estimates at points other than those 
sampled grid points. 

We assume that the antenna temperature measurement at any instrument 
coordinate location x is a linear combination of brightness tempera- 
tures located in the neighborhood of x pi us a random noise element 
denoted by e . The antenna pattern designated by A is assumed to be 
an explicit 2*2 matrix function of x to account for the possibility 
of varying geometries and polarizations as the antenna is scanned across 
the field of view. Equation (1) expresses this rel ationship: 


T A (x ) = //d 2 x 1 A-(x-x‘ ) T b (x"' ) + efx) . (1) 


The instrument coordinate system is defined as follows: 

Xj = direction along the scan with integer values representing the 
location of each sampled data point. 

x^ = direction along the ground track with integer values represent- 
ing the scan line number. 

Unfortunately, the antenna pattern is non -stationary in this coordinate 
system (it rotates with and may undergo compression/expansion as 

a function of depending upon S/C latitude). Antenna pattern 
measurements are usually performed in a space which is simply related to 
"Solid Angle Space", £1 . We shall strongly suggest K -space for antenna 
patternrepresentations in sections that follow. It is usual, however, to 
represent the antenna response function, G , in solid angle space, as 
fol 1 ows : 
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T a (s q ) = / fi(s Ql s) Tg(p ,s) dfi (2) 

where s is a unit vector from the antenna in the direction of the solid 
angle dfl . The region of integration E includes all directions s 
which intercept the portion of the earth visible from the satellite if 
the very small contributions from the cosmic background radiation and 
possible contributions from the spacecraft itself are neglected. In (2), 
G(Sq,s) is the antenna gain 2*2 matrix function (normalized so that 
j ^ trace 153 = 1) and Tg(p,s) is the brightness temperature 

polarizatic s propagating toward the antenna. Tg depends on both the 
region being viewed, which is indicated by the position vector p from 
the center of the earth to the point on the earth's surface, and the 

A 

direction of propagation -s to the satellite. 
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Figure A.l Geometry for Antenna Temperature Oeterminaf ion 
T d corresponds to a two-vector of Horizontal and Vertical components of 

D 

polarization, T * . 
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The antenna pattern components of the 2 X 2 tensors consist of 

S HH * G HV * G VH * and G VV As a first P roblem » ** assume G HV = G yH 
= 0 everywhere so that Horizontal and Vertical components do not mix. 
The generalization to other mixing cases is not difficult. 


A. 2 K-Space Representation 


Our antenna representation is best expressed in k-space (pointing 
vector space). We consider incoming radiation at fixed frequency and 
hence wave vector magnitude (pointing vector in units of inverse antenna 
radius) 


K 2 * * 2 * t 2 * * 2 ■ [f2]‘ 


This is equivalent to having chosen X in units of the antenna radius. 

k is chosen along the “look" direction of the antenna. This (k , k ) 
z x y 

space is important because in this space, the antenna electric field 
response function is the Fourier transform of the antenna aperture 


E = ^dx ( unit source ( k x ,k y)) A ( x ) e1 *' X 

A 

and T^(Sq) » the power, or antenna temperature which is the expected 
square value of electric field <E 2 > is the Fourier transform of the 
antenna aperture function A(x) convolved onto itself. Likewise, if we 
choose R-space as the spatial coordinates to represent the antenna point 
spread function, then the corresponding antenna frequency response will 
be the multiplication of A(w) with itself. Since A(<3) has a sharp 
cut-off corresponding to the finite antenna at jw radius, then so 
will the convolution A(w) x A(w) have a sharp cut-off at jS| * 2« c 

The antenna pattern is presented in coordinates (k , k ) from 

x y 


which 
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k 


z 



and sin 9 = 



+ k 

Tf 


2 

y 


l\ 


can be determined. Also, S = + £ ; S 

0 

axis of the antenna). Inversely; given 



lies along Z (look 
we may construct 


£ = K S n 

0 0 


and 


q = (E - U - <E - EJ 


which ies in a plane orthogonal to S Q permitting the computation of 

k-coordinates for antenna pattern computations. Now, let a be a 

downward directed unit vector along the LAMMR rotation axis. (See 

S * a 

Figure A.2). Then - — = x is a unit vector directed in the negative 

sino 3 

0 

scan direction, z = s Q and s q * x = y can be computed to determine an 
orthogonal set of unit vectors. From this, we can compute 


q x = q • x and q y = q * y 

to give us the coordinates of the antenna pattern corresponding the (S’ , 
S) ; i.e., G(S 0 , S) = G q x ,qy where Gq x>Qy is the representation file 
in k-space. 




Figure A.2 IAMWR Geometry 
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A. 3 Solid Angle to Area Integration 


The transformation from solid angle integration to area integration 
utilizes: 


d« 



dA 


a n 

where p points to the area (km ) element from the earth's center, and S 
is the distance to the area element from the antenna (km). Given the 
position vector of the S/C, r $ ^ c , and the radius of the earth r g we 
note that 


then 

where 

and 


S = S/S 


S 





a 

S 



s/c 


r cos 
e 




where $ is the cone angle between p and r s j Q ; this then completes 
the needed transformation. 


A. 4 Area to Instrument Coordinates 

For conversion to instrument coordinates, we also need to compute 
the area generated by a differential change in the scan pixel (dj) and 
In the scan line ID (di) . 

If v Q Is the ground track velocity (km/sec) and Q q is the angle 
at which antenna scan '.otion on the ground J. to v Q 
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then, r = S„ sin <f> where S rt = r , cos 4 - - r 2 sin 2 f" 

oo o o s/c o e s/c < 


i.s., IPI" * |r„ + $J 2 = r 2 and r 


s/c 


, • 5 = r , S cos 4 

s/c o s/c o y o 


By defining a coordinte system (ground coordinates) via 


A ” r _ / _ -S V— A X A 

5 * ■ ■ ' ■ ’ n = and x = n x 5 . 

r s/c v 0 


then tne unit vector, 4 , along the scan direction at any location 
9 is given by: 


4 * -x cos ( 9 -9 ) - n sin (9-9 ) . 

0 0 


The unit vector , i , in the down track direction along the surface 

a * 

of the earth is orthogonal to both p and x and can be expressed as 


X * p 

A A . 

X * pj 


Here p ~ -l fr-r 2 + r n cos (9-0 ) - r x sin (9-9 ) , 
e 0 0 0 0 0 


- fi 2 

then x * p * n *r * r + 5 r cos (0-9 ) which has magnitude 
e 0 0 'o' 


A * 2 - r 2 sin 2 (0-9 ). So 
e 0 0 


z = n 


2 2 

r - r 
e 0 


r 2 -r 2 sin 2 ( 0-0 ) 
e 0 v 0 ' 


1/2 „ r Q cos (9-® 0 ) 

+ ? / 2 2 " 2 

/r -r 51 n 


r -r sin (9-9 ) 
e 0 v o J 


Now, the angle between i and l is obtained by computing i * 4 
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t x 4 is 


2 2 
r - r 
e o 
TTT2— 777- 


1/2 


r‘ ' r n sin {0-9 ) 
e o ' o ; 

- cos (0-0 Q ), - sin (0-0 J 


r Q cos(9-9 ) 


y_2_ 2 771 2 

* r r q i n 


r ”r sin (9-0 ) 
e o v c J 


x x 


. r cos (0-9 ) sin (0-0 J 


rz ? 

*r ~ r 


r n s™ (9-9 J 


- r cos (0-e ) 
. n -° 0 


ITT 1 


e r o sin ( 0 ‘ 9 O ) 


: ^ r e' r o )1/2 cos < 9 - 9 J 


sin 2 (0-9 ) 
e o v o ; 


+ COS (9-0 ) P 

/r l” r Q sin z (0-0 ) 
e o v o ; 


which has magnitude 


cos (9-9 ) 
o 


X 3 


sin 2 (a-0 o ) 


= COS<A , S> 


where <4,4> denotes angle between l and l . 


(1) The element of arc length generated by dj (one pixel) is 


r A0 
o 


where 


A 9 = Ml JL 

180° N 


5 » __1_ 

6 N ’ 


with N equals to the number of scan points per line. 
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(2) The element of length generated by the s/c ground motion during di 
(one line time = one second) is: 


/l-(-^) 2 sin 2 (0-0J At where At = 1. sec. 
o r o 

e 

(3) The area of (di) (dj) is therefore 

/ r " 1 

d 2 A = v Q A- {^r ) * sin 2 ( e - 9 Q ) At r Q A0 cos <M> di dj 


Finally, 


d 2 A = v AtA9 r cos (9-9 ) di dj 
o o o 


A .5 Transformation From Solid Angie Coordinates to Instrument 
Coordinates: 


The transformation from Solid Angle Coordinates to instrument 
coordinates is done via the area representation as follows: 

t a (S 0 I s I s s> t b < 5 - S) ^ 

A 

= JdA G T 

s 2 B 


Id 2 X G {v n r At A0 cos (9-0)}* ("^) T, 


0 0 


2 ' -B 


where d 2 X = di dj . Hence, A 7 (X -X) 

*o 0 


L±£l 


=- G (S Q ,S) v Q r Q At A0 cos (8-9 o ) 

O 


A-10 




Bi s/\£as i mi Te< hsoum,h: ii. Systems. Is<. 


designates the total »-e1ationship between the antenna pattern in instru- 
ment coordinates (A) and the corresponding pattern in solid angle 
space, G . 


A. 6 The Computation Process 


We begin with the antenna pointing location i Q , j and designate 
X^ , X 2 as arbitrary location variables in instrument coordinates. We 
desire to compute the antenna pattern for all (X^ , X^) when the antenna 
is pointing at location (i , j ) . Then 


T A ^o-Jq) ■ K dX 2 A i ,j < f o- X l V X 2> T B<V X 2> 


Elements to be computed are 

v , r , At , A6 , cos (0-9 ) , 9 , S , p , S 2 and, G(S n ,S) 

0 0 0 0 0 

(a) v Q Computation 

v Q is the ground velocity of the satellite subtrack which depends 
upon earth rotation. It is latitude dependent but can be assumed constant 
for APC coefficient computations. 

(b) r Q Computation 

r Q is a nominally constant radius of the LAMMR ground trace scanned 
circle. It depends upon s/c distance from the earth's center, and the 
earth's radius at subtract location. It also depends upon the instrument 
pointing angle from Nadir. 

r = sin <i> 

00 0 


where 
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S. - r e/( . cos <t> 
o s/ c o 


F- 


r e r s/c 51 " *o 


1 2 

This is computed by looking at S Q = p q - r^ c and p*p = r g - S Q 

r , + 2 r , S cos * . Here £ is the instrument cone angle, 
s/c s/c ooo 3 


{ c) At Computation 


At is the time between adjacent line samples - nominally it is 
At = 1 sec. 


(d) a9 Computation 

A6 is the separation in radians between adjacent pixel samples 



swath 
wi dth 


where H Kusber of samples swath. Nominally, N = 256 or 512 and 
is dependent c.r- che frequency . 

(e) Q Co«puta2i ! '*! 

0 is the angle at which the sweep is orthogonal to the ground 
velocity vector. Nominally, 9 q depends on geometry and on conventions. 
It also depends upon t r.t effects of earth rotiiions. In the absence of 
earth rotation we may establish a geometry for which 0=0 corresponds 
to forward looking; • .e., ( 0 Q = C) . The effects of earth rotation can 
be computed using the transfc- mation equation between inertial and body 
coordinates 


d 

'Xi 


d 

TTt 


- v 
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Hence, v g = v^ - w g X r $ ^ c . The ground observed s/c velocity varies in 
both magnitude and direction dependent upon latitutde: the ground velo- 
city magnitude is 


-V . -v B <^-> 


s/c 


and the direction relative to Vj is 


sin 9 = 


V V B 


0 V B V I 


V I * < w e * r s/c } 
V B V I 


This computation gives 


"etWV - ' r s/c <W 


V B V I 


when 3 0 (circular orbit), then 


r , v 
s/c z. 

Sin “ * r-~ T I " t 

0 NT Hi 

which is latitude dependent. 

(f) 9 and cos ( 9 - 8 Q ) Computation 

9 is the pixel related angle associated with j : 
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0 • % + 40 (j . |) . 

N 

At scene center j - -j and 0 = 9 0 • 

(g) $ and p Computations 

S and p are instantaneous pointings to pixels from the s/c and 
earth center, respectively. They are related through 


r s/c 


+ 5 = p 


We must compute p from our knowledge of i , j Q , and X 9 . We 
also assume initial conditions relative to the earth so that the position 
p corresponds to indices i , j . Then, the change from i Q to X^ 
and j Q to corresponds to a rotation about the ground x coor- 
dinate axis by an angle 


v At 
o 


r 

e 


(x 2 -j 0 ) 


6 . 
J 


followed by a rotation about the new i' axis by angle 

49 < x r' 0 ) 4 • 

The transformation is: 

p . 

o 


cos <$- - sin A. 


sin 5. cos <5. 


1 0 


0 cos <$ .; - sin 5- 

v J 


0 + sin 5 cos 5 . 

J J 
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Note: The order of these rotations are as specified since a reverse 

order would affect a rotation about the x' axis which would not result 

in a translation in the v direction for <5. . 

J 

It is appropriate to expand the 5 . rotations for small angles 

J 

while retaining finite angle (approximations for <5^. . We then note that 
if 


P = T p 


then 


cos 

5 . 

- 

COS 

5 . 
J 

sin 

6 . 
1 


+ sin 

5. sin <5 ^ 

sin 

6. 

1 


COS 

5 . 
J 

cos 

6 . 
l 


- sin 

6 . cos 5 ■ 
J l 

0 


+ 

sin 

6 . 
J 




cos 

5 . 
J 

cos 

5. 

1 

- 

sin 

5. 

l 


+ 

<V 

sin S. 

) 

sin 

5. 

1 


COS 

6. 

1 


- 

<v 

cos 5 . 

( 

l 

, 0 


+ 

< 5 j 

) 




1 

) 


So 


p = p cos («.) - p sin («.) + p (« .) sin (5 ) 

*0 y o 0 J 

P - P x Sin (5.) + P cos («.) - P 2 («,) cos (« i ) 


<\ 3 + P v ( 5 i) + P 2 • 

*o J o 

Further approximations for small 5- leads to the quadratic 
approximation 
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p - P = - P (- 5 -) - P ( 5 .)+p ( 5 . 6 . 

x x „ x „ v 2' y rt v v z rt J i 
0 0 J 0 0 


« z 


' °y 0 * V‘i> - *y„ ( 4) * »z„ «V 


° 2 - V 


♦\ («,> 


The relative error of neglecting quadratic terms is 


«WK [(-£) , 5j] . 


If we are concerned in the APC process only with pixels that are five 
(5) elements removed from the origin, then the relative displacement 
errors will be 


c v At 

max [ | Ae , 5 -2— ] 
c e 

~ [ .025 , .0053 ] . 


This is a negligible fraction of a pixel displacement. Hence, a linear 
approximation is completely justified. Then 
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Computing p from p q given ^ and 5^ permits us to compute 


5 * p - r 


s/c 


(h) $ Q Computation 


(i) S Computation 


(j) [s , s] Computation 


5 o - 5 o - r s/c 


S z = S • 5 


S o " 


5 o 

Fo\ 


s = 


(k) G (S Q , $) Computation 


Suppose we have a file consisting of the antenna calibration data 
it would consist of a two-dimensional array 


Vy 


a ^ 

Given S Q , S , we could construct 

R = K S R » K §, 
0 0 


where 


with 


u , 2" D 
K * y - . 

D * diameter of antenna aperture, and 
A * wavelength of a particular frequency. 
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Then, q = {(R-R q ) - (E-R q ) • S Q S Q } is an antenna related (negative) 

A 

pointing vector projected in a plane orthogonal to S Q . Then, 

S. * C = X sin <fr 
o o 


and 


y = s Q - x 


permits us to compute q = q • X and q = q • y from which the file 

* A A 

lookup in the Gq ,q tables permits us to determine G (S ,S) . 

x y o 


Other useful equations are: 


where 


Also, 


But , 


p , p = Q x p 

o o 


3 = x 5 j + l 5 i 


S - S n a p - p 

O 0 


(S + flxp ) S 
v 0 o' 0 

_ 


and 


A 

s 
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So 


Now 



n * p n (s 2 -s 2 ) s n 

o' O' 0 
T— - $ (S + S Q ) ‘ 


s 2 = S* + ( 0 *pJ 2 + 2 l 

0 ' o' 0 


K) • 


Hence, 


But , 


So 


S - S. 


(G*p q ) s (Sxo q ) 


S o S o 



r 


s/c 


( n *P 0 ) 



r s/c 


(Z*e Q ) 



A. 7 A Preliminary Model of an Ideal Antenna in Instrument Coordinates 

All of the analyses of the previous sections have been oriented 
toward the derivation of an antenna pattern as seen in instrument coor- 
dinates, the coordinate system in which data samples are uniformly dis- 
tributed. In this context, the aperture function of the antenna as pre- 
sented in k-space, if given, would permit one to compute the desired 
response of the antenna in brightness temperatures in the instrument 
coordinates. To accomplish this, the aperture function would be 
convolved onto itself, resulting in the Fourier transform of the antenna 
response function as projected in k-space. From there, the computation 
to instrument coordinates is as delineated in Sections A.l through A. 6. 
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In this section, we consider an ideal antenna of circular aperture, 
having a uniform aperture distribution up to the edge the antenna. The 
antenna diameter is assumed to be four meters. We can compute the elec- 
tric field response of this antenna by considering a unit-point source 
emanating energy towards the antenna aperture with a pointing vector 
making an angle ♦ with respect to the bore-sight direction. 



Figure A. 5 Circular Aperture 
For a unit source, the electric field response is: 


Er - 2 / rdr / d0 e 1 ^‘ r . 

O 0 

- 2 / rdr / de e 1 | k | sin * r cos 9 
0 0 

2J. (k) 

. a £ 
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where 


a = area of the dish = * R' 


and 


K * 2* sin * y 


The power antenna response function is 


A = bE- = b a 2 
Ic k 


2J x (k ) 


1 2 


and must be normalize such that 


/ d a Ag » l 


Hence, 


2 Ji (k) 2 

a 2 / d n(- 1 v -) 2 


we note the mathematical identity: 


; -n » 


.iz cos 9 


J (z) - V / e“ " cos (nfl) d9 


So 


Eg 31 (2)" / rdr J Q (k r sin ♦) 


Let 


kjj =» k sin ♦ , 
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then 


E t = 2. / rdr J Q (k 11 


r) 


= 2 




CdS J Q (O 


we also note: 


f t v J y ^ (t) dt = z J y (z) for Real v > 0 . 


Hence 


z 


/ t J 


0 


( t) dt = z ^ 


(z) • 


So 


E 



Rk n ) J 1 (Rk n ) 


2 * 


R ' 


J 1 (Rk n ) 

(Rk^J ' 


A two-dimensional file was create^ in k-space to represent this 
theoretically ideal power pattern and is the input for further 
transformations to instrument coordinates. It must be emphasized here 
that since this antenna response function has a finite band-limited 
Fourier transform, that an exact representation of this function could be 
obtained by a finite sampling in accordance with Shanon's criteria for 
band-limited signals. This sampling requirement corresponds to a 
sampling of the antenna pattern with spacing <5k equal to one point 
every * divided by two sampling intervals. With this coarse sampling 
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requirement, the antenna function is exactly represented provided one 
does a sine function interpolation between the established values. Since 
the sine function interpolator is computationally inefficient, it was 
decided to sample the antenna function over a finer interval so that lin- 
ear interpolation between points would be adequate. 

In one dimension, the representation of a uniformly distributed band 
limited amplitude spectra with w c = 2 gives the point spread function 
sinc(k) . If we sample this function at uniform intervals k = i Ak , and 
interpolate with a linear Interpolator, we will generate an error which 
is easily viewed in the Fourier domain. The sampled function has an 
amplitude spectra which is the original "box car" spectra but repeated 
at uniform intervals of Aw = . The linear interpolator convolves 

with the comb function of the sampled data to yield a piecewise linear 
approximation to sinc(k) . This convolution in the spatial domain is a 
product function in the frequency domain. The linear interpolator 
Fourier transform function is: 


sine 


2 wAk 

~T 


The choice of Ak in the use of a linear interpolator is determined by 
the maximum transfer error which may be generated. There are two compet- 
ing sources of error generated by a finite Ak . 


• Modulation at the cutoff frequency w c = 2 of the antenna 

• Modulation at the first repeated zone corresponding to 

2 * 

w 3 aJF * w c * 

This is Illustrated in Figure A. 6. The maximum error occurs at w c = 2 
and has magnitude .08 when Ak * .5 . The mean RMS error can be shown to 
be max error » .03b when Ak = .5 . As a function of Ak , the maxi- 
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Figure A. 6 Representation Error of the Antenna Pattern 
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1 2 

mum error Is E_ (Ak) and the mean error is: 

jllaX J 

^mean 3 3/5 * 

For initial investigations, a mean representation error of 3-1/2% was 
accepted and Ak * .5 was chosen. A mean error of less than 1% could 
have been achieved if Ak = .25 was chosen. However, this would have 
increased the size of the required data set by a factor of 4 without mak- 
ing any qualitative difference in the model and/or analysis. In future 
work, a cubic spline interpolator will be investigated. It is believed 
that the cubic spline interpolator will significantly reduce the repre- 
sentation error any may even permit larger values of Ak to be chosen. 

A. 8 Earth Rotation Effects 


There are three major effects as shown in Figure A. 7, caused by the 
rotation of the earth. First, the spacecraft groundtrack velocity varies 
as a result of earth rotation. Near the poles, the spacecraft may be 
moving in the same direction as the earth rotation, thereby causing an 
effective decrease and subsequent crowding of the scan line circles. At 
the equator, the earth velocity adds in a direction nearly orthogonal to 
the spacecraft velocity direction, and a skew of the scan-traced circles 
results thereby. Fortunately, because the scan traces are indeed circles 
of constant radius on the surface of the earth, the local geometry is 
affected in a minor way. In fact, the ground track trace becomes skewed 
in relationship to the instrument axis pointing in the direction of the 
spacecraft velocity, so that the referenced scanner angle no longer is 
directed along the ground track, but forms an angle relative to it, which 
varies with latitude. 

Finally, uecause of this latitude-dependent rotation of the digital 
grid along the scan-line circle and the potential crowding of scan lines 
as a function also of latitude, the set of coefficients that must be 
applied to the data may be latitude-dependent. Fortunately, all of these 
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effects are taken into consideration in the transformation from k- space 
to instrument coordinates through the variable 0 q and its dependence on 
latitude. The maximum angle that 0 q attains near the eauator is approx- 
imately 4°. 

A. 9 Undersampling 

Given the geometric parameters of the LAMMR instrument and the 
antenna diameter, it is possible to compute the required sampling dis- 
tance in both along scan (<5x) and across scan (Sy) directions. For 
an ideal antenna without consideration of signal-to-noise. Table A.l is 
a table of ideal sampling rates for each of the freauency/wave-length of 
the LAMMR design. It is to be noted that the ideal sampling rates are 
significantly higher than those rates planned for the mission. This dis- 
crepancy can in part be reconciled because the real antenna is nowhere 
near ideal. There are a multitude of affects which significantly deteri- 
orate the information content of the received signal. These are: 

• Effects of antenna tapering, 

• Signal-to-noise ratio, 

• Blur Introduced by scanner motion, and 

• Finite horn dimension and placement. 

Except for the effects associated with the finite dimensions and 
placement of the individual channel horns, each of the aforementioned 
degradations have been taken into consideration in computing sampling 
rates for each of the channels that are commensurate with the signal-to- 
noise anticipated of the instrument. It is to be noted that the sampling 
distances computed in Table A.l are based upon the ability to enhance 
signals that are not dominated by noise, thereby compensating for the 
sensor modulation transfer function. We must retain information up to 
that critical frequency if one intends to digitally compensate for sensor 
degradation effects. It is to be noted that even in the presence of 
tapering and sensor motion and noise, that the planned sampling rate of 
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Sq = 1020.4556 km 
D 3 400 cm 


o 



r Q 3 703.726 km 


Ax 



1 .276A (km) 


4y 


AX 

sin 49.94° 


NOISE CONSIDERATIONS 
with TAPER and BLUR 


X 

cm 

' f GHz 

ALONG SCAN 
AX 

CROSS SCAN 
Ay 

ALONG SCAN 

AX 

CROSS SCAN 

Ay i 

.82 

36.5 

1.05 km 

1.37 km 

1.45 

1.71 

1 .4 

21.3 

1.78 

2.33 

2.46 

2.91 | 

1.6 

18.7 

2.04 

2.67 

2.81 

1 

3.34 

2.82 

10.65 

3.6 

4.7 

4.96 

i 

5.87 ! 

i 

4.55 

5.6 

5.8 

7.58 

8.0 

! 

9.48 

5.88 

5.1 

7.5 

9.8 

10.34 

12.25 S 

i 

, 

6.98 

! 

.4.3 

8.3 

11.6 

1 

12.27 

14.5 I 

i 




1 

i 



I 


at 512 samples ax 3 3.6 km 
at 256 samples ax = 7.2 km 
A = 7 km 


Table A.l 


TABLE OF IDEAL SAMPLING RATES 
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512 samples at the three highest frequences and 256 samples at the four 
lowest frequencies will generate undersampling conditions in the 
x-direction for all but the three lowest frequencies, and the same is 
true for sampling in the y-di recti on. 

Further modeling of course, could include other degradation effects 
to the point where this undersampling would be inconsequential. Nonethe- 
less, it is the conclusion of this section that aliasing artifacts gener- 
ated by undersampled conditions must be taken into consideration in the 
design of an antenna pattern correction process, and that such effects 
may significantly hinder our ability to adequately enhance the data with 
context-free matrix convolution processes. Analysis of this sort is 
often utilized to justify placing a low-passed filter prior to digitiza- 
tion, so as to avoid aliasing artifacts. 

We do not recommend a purposeful degradation of the analog signal by 
such a low-passed filter. The only grounds for placing such a filter in 
the system is the ability to perfectly reconstruct the analog signal 
coming after the filter. But, this degraded signal is not the scene. It 
can be shown that a linear interpolator matched to the expected power 
spectra of the original signal (but one which is not necessarily a sine 
function interpolator) will give an interpreted continuous image having a 
lower least-squares residual error than the or.e generated by applying the 
low-passed filter to the original scene. Hence, one can conclude the 
aliasing is not an evil that should be avoided at all costs, but is a 
sign that alternative interpretations of the signal are possible and that 
a priori knowledge of scene content can be utilized to select the best 
among all alternatives. 

A. 10 Antenna Pattern Corre- ^n Philosophy 


In this section the following elements of antenna pattern correction 
are presented: (Figure A. 8) 
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CONCEPTS : 


<> USED TO DESCRIBE A PRIORI ENSEMBLE AVERAGE 
ALTERNATIVE MEASURES OF APC PERFORMANCE 

(1) LEAST SQUARES 

L = <(T b (x) - T b (x)) 2 > 

WHERE T b IS THE BRIGHTNESS TEMPERATURE 
ESTIMATE FROM ANTENNA TEMPERATURE MEASURE- 
MENTS 

(2) , MINIMUM VARIANCE 

5T b = T b (x) - <Tg( x ) > 

L = <(5T b (x)-5T b (x) ) 2 > 

<T b (x)> IMPLIES AN a PRIORI MODEL 

(1) ERGODIC <T b (x)> = LOCAL MEAN 

(2) <Tg(x)> = MODELED 4 PRIORI 

(3) <T B (x)> = STATISTICAL ENSEMBLE AVERAGE FROM 
REAL DATA 

CONCERNS 

(1) TRUNCATION OF MATRIX 

(2) TAILS OF THE POINT SPREAD FUNCTION DISTRIBUTION 

(3) Gibbs phenomena near edges 

(4) alternative statistical assumptions required 

(5) PRESERVANCE OF A CONTINUUM INTERPRETIVE THEORY 

Figure A. 8 APC Philosophy and Concerns 
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• Estimation parameters, 

• Cost functions and alternative philosophies, and 

• Optimization and matrix equations. 

These items form the basis of all the work to follow. The incorpor- 
ation of data interpretive techniques forms the basis of both context- 
free and context-sensitive processing philosophies for antenna pattern 
corrections. 

Estimation Parameters 


In the usual formulation of estimation theory, an estimate of 
Tg(x) , the brightness temperature, is desired. However, for LAMMR pur- 
poses, a more general estimate of the brightness temperature, as seen 
through the "eye" of another instrument (real estate matched Tg(x)) , is 
desired. We denote the antenna pattern of this other instrument as W , 
and the resulting measurement as Tg . The following equation represents 
this measurement: 


Tg(x) = i7d 2 x' Wj(x-x' ) Tg (x 1 ) . 

We desire an estimate of this windowed brightness temperature as a 
linear combination of the antenna temperatures that were measured on the 
integer grid in the x space. The following equation expresses this 
relationship: 


V 5 > - l M x,i V’> 


At this point, we shall assume the existence of some a priori estimate of 
brightness temperatures (two polarizations) over the field of view . Th i s 
assumption is made for three reasons: 
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• The usual formulation of minimum variance estimation filters gen 
erally account for biases separately, 

• A priori knowledge of water brightness temperatures and/or land 
brightness temperatures can be accounted for in this way leading 
to a procedure for removing the influences of land in the near 
ocean vicinity, and 

• The removal of a first-order guess does not limit the applicabil 
ity of the results. 

The following two equations represent the corresponding fluctuations in 
the windowed brightness temperature and its corresponding estimate: 

1T g (x) = //d 2 x‘ Wj(x-x') (T b (x') - <T b (x')>) 

* l M -X,f < T A< ( > - <T A< 1 ») • 
i 

Note that this variance approach introduces the a priori model repre- 
sented by <Tg(x)> . 

There are alternative models that may be impled by <Tg(x)> and 
which stem from additional concepts and/or assumptions used in defining 
the cost function. Three such models are: 

• Use of the Ergodic hypothesis to set 

<Tg(x)> = Local Mean 

• A priori model based upon knowledge of water and land distribu- 
tions on the earth and perhaps a latitude-dependent brightness 
temperature model for both land and water areas. 
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• Statistically acquired brightness temperature ensemble averages 
(stored in a world brightness map) used as an a priori model. 

Other parametric approaches involving pattern recognition techniques also 
exist but are not considered to be within the scope of this section. 
Figure A. 8 summarizes the APC philosophies and concerns with regard to 
estimation concept fo-mil ations and alternative cost functions. 

Cost Functions/Alternati ve Philosophies 

The optimum set of estimation coefficients M should minimize the 
least squares error in the unmodeled fluctuations. We, therefore, 
define: 


x 2 (x) • < &t b u) - «t b (;j) 2 > 

where x 2 is the variance at x and <> represents an ensemble aveage 
over all statistical realizations. 

The cost functions introduced by Backus and Gilbert 2 do not involve 
ensemble averages; they recognize that Tg(x) corresponds to a specific 
point spread function associated with the particular linear combination 
of the antenna point spread functions which generate T^(i) . In partic- 
ular: 


t B (x) = //d 2 x'[jM-jA (T-x 1 )] Tg(x') 
i 

giving the expression in brackets, [], as the actual point spread func- 
tion. Their cost function is: 


2 Backus, G. and F. Gilbert, "Uniqueness in the inversion of inaccurate 
gross earth data", Phil. Trans. Roy. Soc. London, Vol. A26G, pp. 123-192, 

1970. 
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Q x 55 //d 2 x' (Im- jA i (T-x ' ) - Fj(x-x')} 2 Jj(x-x') 
i 

where F is a desired point spread function centered at x and J 2 is 
an associated “penalty function". 

x can be related to Q x by expressing it as: 

x 2 (i) - ////dW;-’ 

(Im- j A i(T-x' ) - Wj(x-x r ) } 

i 

*{IM- jAi(T-x‘ ) -W x (x-x' ' )} 

< 5 Tg(x‘ ) 5 T b (x*')> . 

By assuming <<3Tg(x') 5Tg(x ,l )> = s 2 ^x'-x 1 ') we obtain the equation 
for Q x with a constant penalty function. A penalty function j- (x-x 1 ) 
can be introduced by defining “penalized" "J -windowed" brightness 
variables <5Vg and <5Vg such that if 

x 2 (x) = <(6i? B .59 B ) 2 > 

is defined, then one obtains an expression similar to Q x . 

X 2 (x) = ////d 2 x'd 2 x‘ ' J-(x-x') J-(x-x") 

(Im- ^A-(T- x' ) - W-(x-x‘ )} 

*(lM- jA^T-x' ' ) - W-(x-x‘ ' ) } <6 T b (x* ) 5Tg(x")> . 
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This equation generalizes the Backus-Gilbert method so as to intro- 
duce a priori brightness temperature fluctuation statistics into the cost 
function and vice versa, introduce "penalty” functions into the ensemble 
average statistical approach. 

The introduction of a penalty function is important if one wishes to 
penalize point spread functions having a long "tail" of side lobes. In 
the analysis to follow, we assume no penalty function but recognize the 
tool as another element of APC philosophy that may be relevant near 
land/water boundaries. 


Optimization and Matrix Equations 

2 

Taking a functional variation of x with respect to the estimation 

coefficients M- ?l , leads to the following system of normal equations: 
x , i 

y M- — ( y ** , f + o s,, , ) = r _ •» . . 

s X,1 Vt l ,1 11* 1 X , 1 ' ' 


where 



- ////d 2 i d 2 J" A-O-ii') <«T b (5'J «T b (5")> A-, ) 


and 


r_ . 

x,i 


I I 


J/J/dW-x’ 


Wj(x-x') <«T b (x‘) $T b (x")> A 1 . ll (i"-i") . 


If we now assume that the brightness temperature fluctuations are 
spatially independent, then these equations can be simplified by utiliz- 
ing this statistical independence as: 


<«T b (x' ) «T b (x")> * S 2 a(x I -5 ,t ) . 
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The resulting expressions for Z and r are given as follows: 

= S 2 //d 2 J- A f (f-x"') A-, , (f ' ' -x ' ) 

r x,T" * 52 M* 2 *' W x( x ’-*') A f (f 

Here the £ matrix represents the auto-correlation function of the 
"evolving" antenna pattern in x space; it is an overlap of the antenna 
patterns centered at i and i 1 ' respectively, r , on the other hand, 
represents the antenna pattern as blurred by the window function W . It 
is clear that these integrals can be computed given the analytic forms of 
A and W . If we further decide to select any finite set of sub-grid 
points, x , and truncate the grids i.i' 1 , then the normal equation 
becomes a finite matrix equation which is easily inverted as shown below: 

M = r (Wi )’ 1 . 

It is also noted that this solution minimizes 0 for J 2 = 1 as 
previously demonstrated. 

A. 11 Detailed Approaches to Four Context-sensitive Processes 
• Latitude Dependent Model 


Given world surface map of land flags, ice flags, water flags, etc., 
a priori brightness temperatures will be assumed for each category which 
has both latitude, <£ , and seasonal , t , dependency. For category 
(i ) , we might assume: 

T B (i) * A^ + 8'^ cos 4> + C 0) sin (at+S) . 

Values of A , B , C , a and S can be determined statistically using 
SMMR data. Other terms could be included in both latitude and seasonal 
dependency if the SMMR data can support their need. This brightness 
temperature model would be used to compute a priori antenna tempera- 

A-36 




tures over some fine grid of points which would be stored for each orbit 

in a geographically-referenced world data base for access in real time 

operations. Actually, only sensitivity coefficients y. , to the model 

m 1 

brightness temperatures T„ need be stored since the modeled antenna 
temperature at each point can be computed as <T^> * I y. T g ' ‘ over the 
set of earth feature flag types (i) . 1 

• Parametric Land/Water Model 


In this parametric model, the values of Tg^ are inferred from 

the actual antenna measurements over some local array. The world surface 

map is used as above to compute sensitivity coefficients y^ for each 

point in the local array ( a priori ). Then, we compare the computation 
r { i ) 

estimate <T A > = l y^ Tg v ' to the actual array of antenna temperature 
measurements. The cost function is: 

l [T a (*,m) - <T a (t,m)>] 2 
t ,m 

- I [t a (M) - l ^(M) T B (l) J 2 

4 ,m i 

Taking a variation of this cost with respect to Tg^ leads to the 
normal equations 

A 1 ' i 'B “i ' 
for each category (i) in the region. Here 

B i , « l T a (*,m) y ^ i (*,m) 

* ,m 

and 


- l Y i i (*,m) y i (t,m) 
^ >111 
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and 


is computed as 



A -1 

A i i 1 


Bi- 


Note that typically, the set fi} will consist of one species; occasion- 
ally near water/land boundaries it will consist of two species and that 
rarely does it consist of more than two species; e.g., land, water, and 
ice. An alternative process for estimating <Tg^> is to apply a 
tracking filter and apply the above algorithm only over regions of a 
single species and over which the coefficients are constant. 

• Bi-modal Model 


(T, 


( 1 ) 


In the bi-modal model, two species of unknown origin are assumed 


and 


T 0 ^) . An antenna temperature threshold value are chosen 
B ( 1 ). 


to classify each point in the array as belonging to class 1 {T„ v ') or 

(21 ° 
to class 2 (Tg V ') . Bi-1 inear interpolation is used as a transition 

model between the two fields so that the modeled antenna temperature at 

each point can be computed . The choice of brightness temperatures 

and Tg^ as well as the threshold temperature T h are chosen 

to minimize the variance in the antenna temperatures from the model 

antenna temperatures. The antenna temperature can be expressed as: 


* ,m k ,n 


such that 


(i+A, j+m) e {S 1 } and (i+k.j+nje^} 
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where 


A^ is the effective antenna profile 

$j» are the selected pixels classified in states 1 and 2, 
respectively. 

Sets Sj and are chosen to minimize 

I (T. O.j) - T, (1.J)) 2 
t » j 

* World Brightness Map 

A world brightness map consists of mean brightness temperatures over 
the surface of the earth, possibly segmented by month. This data can be 
used as an a priori model of brightness temperatures to derive the asso- 
ciated antenna temperature needed in a context-sensitive APC process. 


In each of the above models, an estimate is made of the antenna tem- 
perature <T^> which is subtracted from the observed field. 


5t a 


t a - 


<t a> 


The fluctuations ST^ in antenna temperatures are then corrected for the 
antenna spatial degradation in concert with a minimum variance and trun- 
cated matrix correction process whose coefficients are the same as those 
derived in context-free APC processes. The resulting APC produces an 
estimate of the brightness temperature variation <$ T B from the a priori 
model and so one can then construct the brightness tempnrature estimate 


“ <V + 6T B • 
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A. 12 Matrix Truncation Approach 


The system of linear equations th?t must be solved for deriving the 
antenna correction coefficient to be applied at a specific scan angle 
location is given by: 


i 

where the prime on I indicates a truncation over some array centered 
around x . There are two assumption that are imposed in the 
computations that were done. The assumptions are: 

(a) <<5 T b (x) 6T b (x’)> = S 2 <5(x-i') 


. (I- ' 

u i,i 


+ 6 - , 


- » 


x,T' 


(b) W-(x-x' ) = «(x - x* ) 

They result in; 

S 2 // d 2 x‘ A^(T-x’) A- ll (T"-x l ) 
* S 2 A v ,{i”-5) 


with solution 


* 3 ? 9 s - 1 

m - r (I+c 2 i) 

It is appropriate at this time to discuss the two assumptions for 
this matrix truncation approach. In assumption (a), some a priori esti- 
mate is needed for the two-point correlation function in the expected 
deviations of brightness temperature from some a priori model. In a 
context-free case, the a priori model may be chosen as the local bright- 
ness temperature mean over the truncated array, while in a context- 
sensitive process, a priori models may be chosen 
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as dl-cussed In the previous section. Nonetheless, a reasonable choice 
of two-point fluctuation statistics must be assigned. In the work which 
has been done, it is assumed that the fluctuations in brightness 
temperatures from the a priori model is totally uncorrelated from one 
point to the next, and has fluctuation signal power given by S 2 . 

In addition to these assumptions, some a priori windowed function, 

W, must be chosen. Under the philosophy of real estate matching, U can 
be chosen to be any point spread function for which a smoothed brightness 
temperature is desired to be estimated. In the context of the LAMMR 
Algorithm Freeze Report, it is desirable to choose the point spread func- 
tion associated with one of the lower frequency channels as the window 
function for U . In the work done in this report we have chosen not to 
follow this route in order to emphasize an alternative philosophy for 
brightness temperature estimates, which corresponds to brightness temper- 
ature point estimates. For this reason, assumption (b) selects for the 
6 function for W so that the resulting computation leads to the best 
least squares estimate of brightness temperatures at the point designated 
by x . It is to be noted that a broad window function for W will, of 
course, lead to much smoother estimates. However, it is our contention 
that the process of data smoothing should be applied after the non-linear 
geophysical parameters are extracted so as to avoid biases which would 
otherwise be incurred. 

Finally, the matrix equation above which has been truncated, will 
lead to a solution for which the sum of the antenna pattern correction 
coefficients may not sum to unity. This defect results from the lack of 
applying the appropriate constraints in the matrix equation. The con- 
straint, however, is easily imposed by tne following argument. 

The derived matrix of coefficients are to be applied to antenna 
temperature fluctuations. Hence 

6Tg 3 M * 6T^ . 
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Thus 


h m <t b> + M x ( t a ’ <T A >) 

* M * T a + <Tg> - M * <T a > . 

By selecting <T g > * <T A > and writing <T A > = ~ II * T A where II is a 
matrix of one's over the truncated array and N is the number of points 
in the array, then 


t b - " ‘ t a * fr* 1 - h V 11 - T A 

from which it is seen that the derived matrix coefficients should be 
modified so that the mean sum defect (1- ^ . M . .) is uniformly added 
to each coefficient, 

A. 13 Fourier Approach With Alia. ing Consideration 

As previously indicated, it is erroneous to apply Fourier theory to 
this problem because the antenna pattern rotates in instrument coordi- 
nates as the line is scanned. Hence, there is no way to represent the 
antenna measurements as a strict convolution with the brightness tempera- 
tures. However, we make use of the fact that the antenna pattern is 
highly localized, and force a convolution-like expression as a conse- 
querce of these assunptions. The assunptions that must be made are: 

• Ignore rotation dependence of the antenna pattern locally 

• Ignore minor scale changes as a function of location 

• Use a local antenna pattern and its corresponding Fourier trans- 
form in instrument coordinates. 

With these three assumptions it is possible to convert the matrix equa- 
tion to a convolution form as designated below, and to apply the Fourier 
transform results which convert convolutions in the spatial domain to 
product functions in the Fourier d" 
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Basic Equations : 

I H(x-i){I(M") ♦ $ 2 - r(x-V) 

1 

Fourier Transform gives: 

° ' (Mi lA(5)A * (i) + ( I )2 T i l ' W < W > A *(")} 

where »- is the 'box car' function with cut-off at the Nyquist fre- 
quency, and 

M- is periodic: 


•W ■ H ; 

Now, normally, each of the frequency components of the Fourier integral 
are independent. This frequency independence would normally be used to 
generate the system of equations whereby the integrand of the above 
expression vanishes for all frequencies. This would lead to the usual 
Weiner filter which optimizes the enhancement in the presence of noise. 
However, such is not the case in the above equations when the digitiza- 
tion rate is less than two points per wavelength of the maximum frequency 
responded to by the antenna. In such an undersampling situation, the 
Fourier integrals must be summed over all Nyquist zones so as to provide 
a set of independent frequencies so that the integrand can be set to 
zero. When this is done, and the resulting equation solved for the 
Fourier transform of the antenna pattern correction coefficient, , 

we get the following solution which incorporates the phenomena of alias- 
ing. 


£lJ(w+2*I)A*(w+2*Z) 

! 

w l (A(w*2*i)A*(w+2*i)+(|) 2 
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If W £ 1, numerator is y A*(w+2tt£) . 

U 

i 

It Is extremely important to note that the dominating terms in the denom- 
inator of the solution with aliasing at higher frequencies might never be 
the signal-to-noise ratio if the antenna response at the Nyquist fre- 
quency dominates the anticipated noise. A computation of the solution 
with aliasing is In general agreement with the power spectra obtained by 
the matrix truncation approach when the digitization rate is less than 
that required by the sampling theorem. In general, the solution with 
aliasing begs for slight over-enhancement over what would be anticipated 
for a Weiner filter where aliasing is not an issue. These results are 
new and are not, to the best of our knowledge, published in the open lit- 
erature. 

A. 14 Some Antenna Pattern Correction Coefficients 


Figure A. 9 shows sets of coefficients that have been derived by 
means of the matrix truncation approach delineating in previous sec- 
tions. It is to be noted that the optimum set of coefficients drops off 
relatively rapidly with increasing distance from the center pixel. 
Although simulated results have not yet performed on images having the 
characteristics anticipated for IAMMR, nonetheless it is anticipated that 
not too much improvement in the image quality will be seen by extending 
the set of matrices to larger arrays. In fact, frequency response curves 
for these antenna patterns have been computed, as shown in Figure A. 10 (a 
case where aliasing dominates over noise), and those response curves are 
not significantly different. One might infer that the 3x3 frequency 
enhancement matrix adequately compensates for the loss of antenna 
response and that further enhancements with larger arrays may not be com- 
putationally worthwhile. 
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Figure A. 9 Derived Sets of Coefficients 
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Figure A. 10 Frequency Response of Antenna Patterns Corrected 
by Different Matrix Sizes 
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A. 15 Various Error Measures 


The following are various measures useful to the evaluation of 
alternative antenna pattern correction approaches: 

. Noise amplification - It is well understood that any attempt to 
enhance higher spatial frequencies will also amplify noise. It 
is Important that one does not amplify those frequencies for 
which noise dominates the signal. Weiner filtering expresses 
this intuitition very nicely in providing a demoninator contain- 
ing the noise power to signal power ratio. In the context of a 
band limited response, the noise amplification factor is given 
by: 


_ 2 

Noise power amplification (x) s [ H- » 

i 

. Model dependent error measures - if a context-free antenna pat- 
tern correction process is applied, then that process will be 
formulated in the context of a nominal geometric model. In prac- 
tice, the spacecraft will not be in the nominal attitude nr:- will 
it be at the nominal altitude. Both attitude and altitude varia- 
tions will cause a change to occur in the relationship between 
the antenna pattern in spherical coordinates and the instrument 
coordinates of the digital data. The previously derived equation 
gives the explicit geometric dependence of the antenna tempera- 
ture measurement in relationship to the antenna pattern, A . 
Therefore, variations In V with respect to attitude and alti- 
tude can be explicitly evaluated and an error expression result- 
ing from such geometric variations can be computed. The sensi- 
tivity of the antenna pattern correction process to variations in 
geometry has been studied and we have concluded that the effects 
are negligible for APC. 
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. Philosophy dependent error variations - alternative cost func- 
tions measuring the differences between what can be accomplished 
with an antenna pattern correction process in comparison to the 
desired point spread function (with or without penalty function) 
have been proposed. The cost function, Q , is a general error 

measure provided that the point spread function, F , and the pen- 
2 

alty function, J , are normalized as: 

Id 2 x' F-(x-x') = 1 

/d 2 x ' J-( x-x‘ ) = 1 . 

Of particular interest in the philosophy dependent error varia- 
tions are changes that may occur in Q as a result of differing 
philosophyies such as different choices of F and different 
choices of a penalty function. Also of interest are the varia- 
tions in Q when F is chosen to be the effective point spread 
function of another radiometer channel. Unfortunately, time did 
not permit us to evaluate alternative philosophies objectively. 

. Matrix truncation errors - the affects of matrix truncation on 
the error generated in matching the desired point spread function 
can easily be determined from the expressions previously 
derived. In general, a simple truncation of coefficients is 
inadequate because the sum over all of the coefficients should 
add up to one in order to preserve the appropriate response at 
low frequencies. Therefore, it is necessary to redistribute the 
truncated amplitudes to interior elements of the sum and it is 
possible to do this in several ways: 

- add it to the center element 

- distribute it proportionately to the coefficients 

- distribute the differences uniformly over all of the elements. 
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As would be expected, different optimization philosophies can 
lead to any of the above. Within our work plan we have concluded 
that the third option above is appropriate for the statistical 
assumptions made in the derivation. 

• Uncorrectable errors - it is virtually impossible to correct com- 
pletely for the blur introduced by the antenna pattern. In the 
context-free situation, nothing can be done about frequencies 
that lie above the effective cut-off frequency of the sensor. 

This is important to understand because the sharp water/1 and 
interface is a discontinuity containing many high frequencies. 

The BTS proposed measures of these errors incorporate aliasing 
phenomena. 

* Context-sensitive improvements - Section A. 11 and 3.1 contains 
our presentation of the incorporation of a priori knowledge in 
brightness temperature estimates. In a context-free scheme of 
antenna pattern corrections, nothing can be done about spatial 
frequencies that are higher than the antenna cut-off frequency. 

In a context-sensitive antenna pattern correction process, the 

a which has been gathered by the antenna is analyzed and 
interpreted in the context of our a priori knowledge of scene 
variations. BTS envisions a context-sensitive processing scheme 
in which the high frequency content of brightness temperatures 
are inferred from the measurements and our knowledge of the real 
world. The improvement in brightness temperatures which can 
result from a context-sensitive process must be determined 
through simulation; this work requires simulation and has not yet 
been completed. 
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Appendix B 

FORTRAN program listings of the APC demonstration software 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


SUL^QUTINE LENS(RAD ILS ) ... 

C 

C********* ******** ** 44 ** 4 ** 44 * 4 * 4 * 44 * 4 * 4 ** 444444 *********** ****** 

C TO GENET: ATE A £4 > 64 AIRY DISK, WH I CH CCRRESPCNCS TC The 

C SQUARE OF THE FIRST-ORDER BESSEL FUNCTION DIVIDED BY ITS 
C ARGUMENT . 

C US IMG THE CONCEPTS QF LINEARITY, THIS IS THE IMPULSE OR POINT 

C~ SPfiEAC FLECTION OF AN - 1 NCChEREN T . S-Y S 7EM CUE TO A CIRCULAR 
C APERTURE WITH A 0-1 REAL FUNCTION HAVING l IN THE AREA CF THE 
C APERTURE AN D 0 ELSEWHERE. 

C Tt-t: IN-TOUT Cl SPt. A T--CCNT-A-1 NS TH^ MAIN- tCee ONLY . 

C AIRY DISK IS STORED AT AlRY.CATA CN CISKXX USING UNIT » It 

C 

C***** ****** ************** 4*44* « 4 * 44444 4 ********* ***** ******* ***** 

C 

DIMENSION AIRYU-4...6A) ... A( 641 

V* 

GR ID = RAO I L 2/S 2 • 

- Xt ~ -PHAOi-LHS * - 

SUM = 0.0 

REAL = l .0/( GR ID*GR ID ) 

c ‘ - - • 

DO 10 J=l.£4 

. U = U+GRID . 

V - RADIUS 
DO 10 1= l , 6 4 

R = SURT(U*U+V*V) 

IF(R.LE.3.QJ GO TO 33 

ANGLE = ~ R - 2."2561SS + 0.124fS6>3./R >' 0 . 000 056*9 ./ C R ) 

+ - O.CC637S*27./( ft* R*R ) + 0 • 00 0743 *8 1 ./ ( R* R* ?* H) 

EL-= .a .7-S.7£E£ * C-CCCCC47/R. 4 Q» Q16S9? *S •/ (R *R 3 
+ H C»000171*27./( R*P. 4R ) - C.CC2495481. / {R*R*R*R ) 

SC = FI * CO S ( ANGLE )/(R**1.5) 

*-ifl*d-r3-F ■— EE4BC44-. — 

SUM = SUM 4 A IR Y ( I » J ) 

_GC TO 10 _ 

S3 BC = C.S - C .S£2450*R*R/9. * C. 2 1 C936 *R *;*R*R/8 1 . 

* - -D*A39*.43'*R.*JR*ft.*S.4<R.*fiy 72S^-* 0 •.0Q44D3* C R* *8 )/ 65o L » 

A IRY ( I • J ) * eC* 8 C* 4 . 

SUM * SUM ♦ AIRY ( I , J) 

10 C&NT-INWE 

C 

W R IT E ( 6 , 24) SUM 

3CAL.E = REAL/2LM “ ‘ * ~ 

RSUM = 0.0 

DC L^.6.4 __ . .. .... _ 

DO 23 1=1,64 

AIRY(I.J) = AIRY(I.J) * SCALE 

R 50M -—ft-SW — - 

22 CONTINUE 
22 CONTINUE 

- jyr ite rc , in ) — 

171 FORMAT ( * PRINT POINT SPREAD FUNCTION 1 (NC=0)») 

. . _RE>DCE. 1.72) IPRXN.T ._ . 

172 FORMAT(Il) 

C 

-HH WTINT^EO-^)- GO- - 

C 

WRITE (£,24) RSUM 

‘24 FCRMATC SUM OF* AIRY" = * »F5"i ' * 

DO 12 J=25 ,39 

.. DO .L2_I=25»29 _ ...... 

A(I) * A IR Y( I , J) 

13 CONTINUE 

WRITE! 6-r-17 3-) — T-A(-I )- 39 H - - - _ . 

173 FORMAT ( 1E( IX ,F4.2 ),// ) 

12 CONTINUE 

*RITE( £, 174 ) 

.174 .FORMAT. (///.). . _ 

C 

c 

d -©G- -IEC- U*-4-r-£-"4 — 

*R ITE ( 12,01) ( AIRY{ I ,J) ,1 = 1 ,£4) 
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^ao — ccmzuiuL e 1 2 * 6 i 


RETURN 
~SHi 


C T!iT - DRGCS * M -?°R .SUBROUTINE .IfNS 
_Q_ _ DIMENSION AIRY( 20 .20 ) 

** ite c e7 rr — 


^ KM AT(* EN TER rat tie — _ 

T > Di SI 

.2 PORMAT(F 4 . 1 ) ‘ 


"CALL LENS CBTADTVS) 

• S.IOP 

END 
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ORIGINAL 

OF POOR QUALITY 


C MAIN PROGRAM : AN TE N A 

C 

£444 4444444444 44 444444444444 4 44 4 44 « «<< * * *44* 44 *4** 4 4 

c 

C TO 3ENER ATc A FILE ANT. DATA CF <4X64 ELENENTS 

C — V-AL-Oes rN--!MSTfn;MEN-T -CCOR^trfA-TEJ-i "THAT I S“f— THE" X 
C SPACING IS 7/4 XILCMETER. PROGRAM TESTLS SHOULD 

C IN CROLR TO GENERATE THE IL'EAL ANTENA PA TTE P AIR 
C J AT A 15 STORED ON DISKXX BY I /C LMT « IP * THE 

C WILL 8E STORED CN DISKXX dY I/O LN I T A 1 1 <. ANT. 

C MENT COORCINaTES.. aHILS A IK Y .DA TA I S I N O- SPACE ► 

C TO THE POINTING DIRECTION. 


4* 4**4444 444 4 44 4 4» 

CONTAINING ANTENNA 
SPACE. THE GRID 
BE EXECUTED FIRST 
'1*4.64). AIRY. 
RESULTING ANT. DATA 
DATA IS IN IN3TRU- 
A PLANE ORTHOGONAL 


i 4*4'4*-»4r4'44*-4*44r4*4r- 


CQMMGN AIRY(<4*o*) 
DIM2NS ION A ( < 4 ) 


AN T( 64 .64) 


-Fuxssu- — - 

FJY=0. 

NO SP 3= 256 

OK ID -« 0~ * 2~S 

SUM « 0.0 

REAL = I .0/1 GR ID4GRID ) 


WH ITcK. 30) 

311 - EORMA.T C/-w* -EJvTER- -FREU-UENOt 

REAL (6.31) GHZ 
31 FORM AT IF6 .2 ) 


IGHZi- L*) - 


00 280 J=1 .64 
READ! 12.31 ) (A IRYt I 
280 CONTINUE" " 

F I Y= — 9 . 

OO— LL -Us-U-EA 

FIYSsFlY + GR ID 
FIX=e . 

S O H fc- I - l . C* 


J) .1 = 1 .64) 


F IX=F IX-GR ID 


"CALL TXF'OfTMTF'J >.FJ Y^.FT* ,F I Y7 CX . C T75FTZTKCSFS .CC^S T > 
CALl. INT3P (QX.QY » A X .A Y .8 X ,8Y ,CX ,CY ) 


-CAU IN-TLNXOJL^ax. AX .A Y » flX ►EX , 


g NTI I .J)^ = 


VALUE 
— AN- f -t-i 


4 CONST 
tUT - 


11 

10' 


CONTINUE 

'CONTINUE" 


*R1TE-1£..2A)„SUM 

SCALE = REAL/SLM 
RSUM = 0.0 

Oe— a-2— > *-Hr-64 

DO 22 1=1*64 

ANTI I . J ) = ANTI I .J ) 4 

RSGM“=~T r 5UM 'ONTrn'j : 

CONT INUE 


SCALE 


23 

_ . 22 .- 

1 11 


CQNT LNLE 

nRITEIE. I ll) 
FORMAT!* PRINT 


1 12 


ANTENNA 

REA»fSrL4 a) I P- r . IN* 

FORMAT III) 

IF! IPR 1NT.EQ .1 ) GO 70 113 


PATTERN i l NC = 1 ) ?*) 


WRITE 16. 24) R SUM 

24 FORMAT.! !_S.U M__QF_ANT _=_ 

DO 12 J=25.3S 
CO 13 I=»l.«4 

13 CONTINUE 

WRITEI6.173) I A{ 1 ) . 1 = 2 5 . 35 ) 

TT3 — TERM ATI T5TF5“.4 ) ./7T 

12 CONTINUE 


£‘^2i/.Vi 
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lU 00 lcC J=* l • 64 

*£..ITE-(_LI +&.I.X- t-Att-T tT *4 - X ,1=1^4) 

81 F0-XATI64dl2.<«) 


18MAR81 1 i « i'3 • 22 - VOL=SCFIOOO. DS*=YC* j£. F CLA NO. FORT 


180 CONTINUE 

c 

»R ITS i C.200) 

20 0 FORMAT (/«* THE V -A XI s 1 » } 

201 FORMAT! l0F8*O 
MR IT5(e, 202 1 

73TZ~ FORMAT (/"T * 'TFe~X^A*>TS~l'» ) 

4RITE(6. 201 > ( AN7( I ,32 ) • I *32 *£+ ) 


STOP 

EN3 


... At 




3U3R0UT1NE 7/FCRMi FJX.FJY ,F I * ,GY . GhZ . NCSP S .CCNST ) 

C. 

C«** ********** **4r* ********** 4 *** ** 4***44**<# ****4*** 4* *4***** ********* 

c 

<• TO -TRanSFO*** -GPGt^O—iNS^RtfWeNT CCCPOi-NATeS DEARTH SUPf-*Cc) TO — 

•- ANTEMNA COORDINATES (A PLANE PARALLEL TC ANTENNA APERTURE). 


C 

C — 

c 

c 


INSTRUMENT' COGRCINATES APE" IN ‘RF I XEL NLMEcR AND SCAN LINE NUM3SR . 
(FJX.FJY) : POINTING DIRECTION 

CFIX.jr.LY.L_ I ANY. NtU.Gt'BGR- POLNTS — 

(QX.QY) : ANTENNA COORDINATES 
GFZ I PRIMARY PAOICMiiTER FREQUENCY 


CONST 


ANTENNA CONSTANT 


C* ********** ** *T* *:**** * * *~4 * * 4**> * ♦ * **4*4* ***************** *** ** ***** ** 


C<<<<<<<<<<<< 


constants 


«>Nxn*\\\Nv 


DIMENSION RSC£3) .0P( 31 . PC( 3 ) .P (3 1 . S (3) *SC (3 ) .US (3) .USC (3) .VK( 31. 
+ VKQ ( 3 ) .DU'" 2 ) • XX (3)»YY(3) .SNAKE (3) .Q<3 ) 


PI = 3.14159365 
REARTH = 6371.0 

■“ OTIME'= IVTT - -• •' 

DThETA = (PI/l£C.)*l5C ./FLO A T(NCSFS) 

TEETAQ_J=- Q..Q. 

VTRACK = 7.0 
ALTITO = 700.0 

RANNT - Q. QC3 
C - 3OCOO0. 

ThET - A i OTHETA^'FTX ~ " ' 

C 

C <.«_«_ ..DEF INE SP ACJLCRAFJ_gJLSIXICM-._FEE-OUE.hGY-. ANDL-P-CI NTlNG DIRECT IQM >. 

C 

P.SCl 1 ) = 0. 

r»S-£ (-£-)- » — 0-. : ~ - 

WSCC3) = -ALTITU-REARTH 

AESRSC = SQRTCRSCC l ) *R SC C 1 )*RSC ( 2 1 *RSC (2 ) +RSC ( 3 J *RSC(3) ) 

ARS C PO = r~ARSTNi rAL:TI TC*REAHTFT«Sr Nr PBICr7TTEAPTH)~PHIC 

PQCl) = 0. 

PQC-21 =_3EARTK*SIN C.ARJlCP Cl 

PO(3) = -REARTF*COS(AR SCFC) 

F = GHZ* 1.0ECS 


C<<<<<<<<<< TRANSFORMATION ) DD > J ? > 2 > > D Q 


ITT^ D1T- ET A * C F nr-T'JXI 

DJ = V TR ACK *CT IME/REAR 7H* (F I Y— F J Y) 


C 

ELl-L = P. QL LQj^PCI. 3-1*0 1 — 

P I 3 ) = PU(3)*PC( 1 )*0I-P0:3)*CJ 
P ( 3) = POC 3)*PC( 3)*OJ 

— c 

AB3UP = SQRTtP (1 )*PC l)+P( 2 ) *F (2 )+F{3 ) *F(3 ) ) 

UP ( 1 ) = P ( I ) /AESOP 

UPTZT = PT"Z LTTE'S OP 

OP ( 3 ) = P ( 3 ) /ABSOP 

.C 

S( 1 ) = PCI )— R SC ( 1 ) 

S ( 2 ) = P ( 2 >— RSC ( 3 ) 

3< - 31 - P -( 3> — RSCC 3) 

C 

S0( i ) = P0( 1 >— R SC( l) 

SCTT J = PO ( Z J-TTSCl 2 ) " ' 

SO (31 = POC 2L-RSCC3J 

~ abss = sort c sc u*s( i ) Vs { 2 ) * s (TT+s IsT »sTf f> 

ABSSO * SORTCSOC 1)*S0( 1)+SC( 3 > * S C <2 ) +SC ( 3 1 *S C (3 ) 1 

-€ 

C 

OSC 1 ) = SCI ) /AESS 

— trsr2)'=~s(2)/ABss 

US ( 3 ) * SC3)/AESS 

usoTTj ”= so c i Ty absso ' 
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GSOt i ) = sot 1 ) /AB3SC 
— usut j i ^.sat^: wales C- 


• AVfeL = C/F 

F K — em *(H » R A NN-T->H»AV feL - 


VKt l ) = FK*tSt 1) 

VKT2T "= FK TVJS C TT 

VKt3J * FK*LS(3) 


VKOl 1 > = FK*CSQi 1 > 
v KQtjl - ric *!i" ci") 


OUMt I J = 
OOMTZ T”= 
OUM( 3 > = 


VKt l J-VKOt 1 ) 
'VKTTT-VKGrZr 


VKt 3>~VK0( 3) 

_SaoaiK_=_ USO t LJ^xlUKLJLX+LSat2J-*0 LM(2 L+LEC C31*CUUF uij_ 
Q ( 1 ) = OUMC l >— SQ0C7K*t. S0< 1) 

Qti) * OUM( 3 )-S0C0TK*uict3) 


SNAKE CTJ =—-J 

SNAKE (2) * -RSCt 2J/A6SRSC 
SNAKE, t 3.) a -R SCt 2 > / ABSRSC 

ANGLE = SIN(PHIO) 


-* X« t - K — H> 6 Q " t 2 - ) » S N AKE( g» - i : S0(3) »SNA *L F2- ) - W AN GL E- 
XX ( 2 ) * (U SO ( 3 ) * SNAKE ( 1 ) — l. SO 1 1 ) A SNAKE ( 2 ) ) / ANGLE 
XX (3 ) = (USOt 1J* SNAKE l 2)— CS0t2) ♦SNAKE(l) WANGLE 


YYtl) = USOt 2>**Xt 31-ISC(3)*>»(2) 
-YYtZ) _=_USQXJL> »XAt lJ- L £Ot.U»Jr JUJU. 
YY ( 3 ) * USOt 1J*XX< 2)-LSO(2)*>X( 1) 


GX- 


QY = Q( 1 >*YYt i )+Qt 2)*YY{ 2T+Qtl) ♦YY(IT" 


~RO ”=' 

CONST = VTP.ACK * RC * DTI ME * 071-ETA * CCStTHETA— THCT AO ) * 

tU SLU.T LP.t l ) + L St J2J-*U Fi 2 } + US LI ) * L F 1 3 1J S * ASSS ) 

♦ * FK * «QRTtFK*FK-VKt l) *VKt 1 J-VK(2) *VK(2) ) 


Q CTi iCki 

* ■ vA W 

END 
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SOfcHG 0-7 H*E— I H- T r P 4 fr V -,-**-. AY T&^r&Y — 

► +♦**#***♦♦♦* ************ ** * * **** 4 * ** ♦♦ ♦♦♦ * 4* ♦♦♦♦ ********* 4**44 *♦♦** 

TO LOCATE 2 NEAREST POINTS TC <C*,QY) IN THE Q SPACE. The POINTS 

.ARE. QM TirE 41.5 SPACING GR1J.. lh£l£ 3 _ FCINIS WILL QE USED TO ... _ 

GENERATE A PLANE FOR IN TEF FOLAT I NG THE VALUE AT ( CX , CY ) . 

-lUX, 0-V->- IMP OP- POINTS — - 

A.B.C ThPEE NEAREST POINTS 

*************** ******** ****************************** *************** 

DIMiIttSJON PX( _ - - 

C 

IGX * CX 4 SIGN(C.S.QX) 

-H3 a ■ F LO A T* I Q->4 

S = QX — FQX 

QX1 = FQX 4 SI GN( C.E • S ) 

(TXZ - ~ F<TX ~ “ ' 

c 

igy, =_-QY— ». -SIGN! JL. E ^OYJ . - 

FQY = FLOAT { IQ T ) 

S = QY - FQY 

OY 4 -- » -R 4Y — 4 --SIGN <- 0 ~»e.e ) — 

QV2 = FQY 
C 

pjcm-cx-r 

PY(l)=QYl 

. PX (.<1=0X1 

PY(2 J=QY2 
PX ( 3 !=QX 2 

PAM-34 = - Q* - l 

PXC4 j=QX2 

_PY£41*QY2 

00 IQ 1= l . 4 

MUCH (_L) _ = _AaS.tP X(-I_J~gXL-f-Afl.S.(P.Y-LLl^Q Y 1 _ . 

10 CONTINUE 
C 

&-1-G — * — G-»-0 — - ■ — — - — — — ______ 

00 11 1= U 4 

IF(MhZCH( 2 ).LE.dXG) GC TC 11 

'HiG=whTCHrrr _ 

MAX= I 

11 CONTINUE __ 

C 

1= l 

TP47-.e Q.HAA ) — I— M-4 — 

AX=P X ( 1 ) 

A Y=P Y( I_) 

IF(I.cQ.MAX) 1=141 

„.BA=PX(.JJ 

8Y=PY( I ) 

1=14 1 

RF-t- I .EG , HA X -)- - I « 4 -4-l 

CX=PX( I) 

CY=PY(J)_ 

RETURN 

_ END _ . 
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' SOeRGUTINE INTLN(0X iGY .'AX ,AY .BXVeY .CX.CY .VALUE) 

c 

C.********* **************** ********* ******* **********•*****-*****♦*******• 

GIVEN 2 POINTS A. 0. S C. LOCATED CLCSE TC Q. ANC GIVEN F( A) . FIB 

— &~~F f€ >» TH-IS- -S L -B W fc-T-IN C iN T€RPC-fc*-Te-S- F-IG> -tf-NEARtY • YAttJcr -= F tQJ • 

F ( A ) * c ( B ) * AND F (C ) ARE CBTAINED FROM AIRY(X.Y). 

*** ************************ ****** *** ***** ****************************** 

. . COMMON. .A IR-Y tt<L*A>4 J 

DIMENSION A( 2 ) ,e<2) »C( 2) .AS( 2) .AC (3) .A£XAC(3) 

C 

t A as i — t-5 t *r-** irrOv^Hrm 00 

INDY=1NT(A3S(- 15. 5— A Y ) / 0 • 5+ 1 -0) 

A ( l)*AX 
A ( 2) = A Y 

A( 2)=AIR Y( INCX • INDY ) 

INOX= INT ( A3S{-ls^5-BX )/Q .5-M -0 ) ~~ ~ 

IN0Y= INT (A3S{— IS *3— B Y ) s 0 • 5+ 1 .01 

e( 2)=BY 

_ at Z) = AIRYt INQX, I NDY ) . _ 

INOX= INT ( ASS (- IS. 5-CX )/<}.*+ 1 .0) 

- _ .INQY=INTCAflSt=r 1E.S-C.Y Lza^«.+ I .0) - - 

C< 1)=CX 
C ( 2)— C Y 

F— A-1.R-Y <- INO-X- rH«t>Y-i 

A8( l )— A( 1 )— B ( 1 > 

ab(zt=ac zy-erzj 

AE(3)=A(3)-B(2 ) 

AC t Li= AC Ll-CLLJ --- 

AC(2)*A( 2)-C(2) 

AC13)*A( 3)-CC2 > 

ABXAC ( 1 l=A0 ( 2)*AC(Z)-AB{2 )*AC(2) 

ABXAC(2)=AB( 2 )*AC( l >-A e< 1)*AC(2> 

AFX A C (“3 J— A BT l JAXCCTT^ AFTZ) *ACL I ) 

C 


VALJE * ( — AE XAC( 1 ) *< Q X—AX )— AEX AC ( 2 ) * < CY— A Y ) )/A£X A C ( 2 ) + A(3> 


RETURN 

END 
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C TEST PROGRAM FOR CCNVOL AND SIG1A A NO GA*NA AND CCE F ANC P5F 

COMMON ANT (t4«(4)i LMAT IX!49.49) •GMATI>(4 < 5) »P SF (04 ,6 4 ) 

_ DIMENSION I A ( 49) > •FRLANTJ, 64 .,64 ) ,F I RANT (64 ,C4> _ 

C<<<<<<<<<< READ ANTENNA PATTERN FF CM I/C » 11 55>555>>>> 

_ JO 280 J= l ,C4 

RE: AO (11*81) ( ANT ( I . J ) . 1= 1 *64 ) 

»i FefrM^TltrAeirZTrfr) ~ * 

280 CONTINUE 

C ..... 

C<<<<<<<<<< ENTER APC PARAMETERS 5:55::55>5 
C 

. WRITE! £.-1.01- - 

10 FORMAT ( * ENTER SIZE (7. 5, OP 2> .*’) 

RE AO (Still I SIZE 

— F &SMA-Tti-t-) - 

WRITE!*. 20) 

20 FORM A T ( 1 ENTER SQUARE OF THE N/ S RAT I C _( F6j»4 )__:•) . _ 

REA015.2IT" SNTSR ‘ ' 

21 FORMAT (F 6*4 1 

* jutec £.. 3 .ii - -- -- • — 

21 FORMAT! • PRINT MATRICES ? !NC=t) : ') 

READ ( 5 *22 ) IPRINT 

. T 5 g -*^ 4 T i 1 f - — — — — *■ - — — — 

■ QJ “UnW ^ I f T I / 

MSIZE = ISIZE4ISIZE 

c _ _ 

CC<CC<<<<<< COMPUTE SXGM^ AN C"OAM»TA- MATRICES* 5 5 5»»»> 

C ^ C TfSMAf T «T7P l 

<T 

IF! IPR INT*EQ .1 ) GO TO EO 

w ft-iTE-t-*-.-*-S 3 

CO 12 J= 1 • M S IZE 

CO 12 I=1.MSIZE __ 

ATTT=EW*TTXnT J) “ ~ 

12 CONT IN US 

WRITE!*, 14) !A( I J *1*1 .MSIZE) 

_.1A FORMAT !■&! 1 X » FE^A-1 +J - 1 

, 12. CGNTINUE 

15 FORMAT!///) 

-6 

50 CALL. GAMMA! ISIZE) 

C 

imp RTNT .ECTT 'TO *T 

CO 22 1=1. MS IZE 

A.L.L1 = GM AT I >t I 1 

22 CONTINUE 

WRITE!€,14) (A! I ) .1=1 .MSIZE) 

<r 

C<<<<<<<<<<<< SOLVE LINEAR EQUATION F CR APC CCEFFICIENTS >>»»»>»> 


WR ITE ( £ * 34 ) 

FORMAT!* THE APC COEFFICIENTS :*) 


'aR IT E ( 6 . 15) 

C 

C«C« < ««<"C~COMPLTE -PCrNT* SPREAO ‘ FOCTT CN T5 55555 5 5> 5 5> > > 
C 

CAU ANTPSF! r?t7g] 

C 

WRITE!*. 41 ) 

-4-1 F O RM AT F/ -i «— P R I N T Pgr t ~ ! N C— t ) — ?-»• ) — 

RE AO ( * . 33 ) IPRINT 

IF! I PR INT. EO . 1) GO TO 42 

WRITE!*. 35) 

.25 FCRMATX»_ THE. POINT gP RFAP EiiNCT-I-CN L* J 

00 22 J= 25 , 39 

WRITE!*. 172) !PSF! I ,J ) ,1 = 25.29) 

1-7- 3 - FORM A-T-(-l-SF-S-w-S-r*“) — 

22 CONTINUE 

wfi IT£ ! *. 15) 
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*R ITE * 6 » 175) 

175 FOSMATC THE Y-AMS : • ,/ > 

rtfflTECf t T74J (PSr(2ZrS) , 3=32 .64 ) ‘ 

174 FORMAT* 10F8.4 ) 

WRITE-Cti L76J. _ 

176 FORM AT (• THE X-A A I S :••/) 

MRITEC6. 174) (PSF(I. 32). 1 = 32.64) 

C 

SUM= 0*0 
00 1 J=l.64 

'CCT2'I=IY64 — - 

SUM=PSF* I.J)+5LM 

___£ _ .CONTINUE 

1 CONTINUE 

MRITg* 6, 177) SV.M 

H7-7-FGRMA-7 < » - T HE SUM O F P SP — - * .F- 7 ^ 0 -. ^ )— 

C 

C<<<<<<<<<<<< COMPUTE 2—0 O FT > > 3> ) 3 3 2 3 >>> 
42 CALL CFT2D* 64.64 ,P«F .FRLANT.FIFAM) 



45 


46 


M R IT E ( £ * 43 ) 

FORMAT*/. • PRIM OFT 7 * NC=l ) :*) 

R E- A M-frr^J-)-- I PQ INT 

IFdPR JNT.EQ .1 ) GO TQ 144 


WRITE r£.“44 I 

FORMAT*/,* THE Y-AXIS (REAL) I',/) 

*R IT 6 ( fj l 74 1 ( FRLAMI L22 ^J±,J-22L. 6.4) 
MR ITE ( 6, 45 ) 

FC RM AT ( / » * THE X-AXIS (PEAL) : * ./ ) 
-MR4T-6- 


MR ITE (6.46) 

FORMAT*/, • THE Y-AXIS IMAGINARY !»,/) 
WRITE* 6. TT4 ) t FTMANTT ZZ, XT. 2^2Ti 64T ~ 
MR ITE* £• 47) 


47 FORMAT*/.* THE X-AXIS IMAGINARY :*,/) 

- WRITE ( 6»L74) L FLM AN-TLL-,32 ) W =3A^4J 

C 

MR ITE * 6, 46) 

46 FeR*4AT</-,-*— THE-RE/Ht OF-T — *%/} 

00 150 J=25.25 

CO 151 1=25.25 

A Ci ) =FRL AN T * I. j ) 4 256.0 
151 CONTINUE 

MR.ITE-t-6._L73-)- LA* I- ) .3=3 5.35-) 

150 CONTINUE 

MR ITE *6, 49) 

49 — F O RMA T* / . »■ Tt»E~* IMAS T NA F' f — *V/~) 

00 52 J=25 .35 

_ 00 53 1=25.35 

A * I ) = F IM ANT * I. 3) *256.0 
53 CONTINUE 

.... W p J T P ( F r 17T ) * A* I) 

52 CUNT INUE 
C 

~C 

144 STOP 

__tNCL 
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oUbRULT I NE SIGM A f ISI Z E ) 

C 


c 

C TO COMPUTE THE £ MATRIX. S REF F ESE NT S THE AljTO-CCFR ELAT I CN 

C FUNCTION CF THE ANTENNA PATTERN T N~ THE X-SPACE.' " TT'I S AN OVERL AP 
C OF THE PATTERNS CENTERED AT ANY 2 FCINTS fcITHIN A TSIZE X ISIZZ 

C JK INDQHi. * HERE. ISI Z£-_= IHE.-E MATRIX.LS THEREFORE. A 

C 49X49. 25X25. OR 9X9 MATRIX. THE CFDEF CF THE ELEMENTS Hr ITHIN E 
C IS DEFINED BY IPG INT AND JPCINT. 

C*** ******* ********* ******** ****** 44**4*** **************************** 

” CD M7TO N _ A NTT6"4'?*E'4T i EMXTTXT^S V4 9TTGMATTTX CA ST .PSF(6 4^,64 T 

DIMENSION IPOINTC 49) . JPOINTC 49) 

MSIZE^IS IZE*I SIZE " 

C 


ISTX= ISIZE/2 
I STY=— ISIZE /2 

DO 1C 1=1. I SIZE 
-U5.G1RT Cl LI js_l.£2Jfc-.l4L 
II = II+l 
CONTINUE 


1 1=0 

otr tz - j= r, rs izte - 

JJ=ISTY+J-l 
.00 li K=-1 aJLS1ZE_- 
1 1=1 14 L 
JPOINTC I I >=J J 


12 CONTINUE 


DO 14 J= 1 .MS IZE 
1X2= IPG INT ( J ) 

lY2=OPaZNT<U-)- 

00 15 1=1. MS IZE 

IX 1= IPCINTC I ) 

— IT t=-jf*tfl NT Mi 

: 

CALL _C0NV0LJ i:.t» IV l_,I >2 ,I_Y2 ,EI1 12 ) 

EM A T I X ( I . J ) = E l 1 1 2 

c 

— 15 CONTINUE 

14 CCNTINUE 


Oc Z'" L P ^r is 
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SUBROUTINE CCNVCLt IXI.IY1.IX2.M2.EIU2) 

C*m* ****** ********** ******** 4 + 4*4* 4*4 4*4 4 4 4 * * * 4 *4i 4> **4^ * 4* * * ** **** *** ** * 

C - TO -COMPUTE - AN ELEMENT rN~THE‘~E' _ MATP’I X. E REPRESENTS THE AUTO- 
C CORRELATION FUNCTICN CF THE ANTENNA PATTERN IN THE X-SPACt. IT IS 

C AN OVERLAP OF TEE ANIENliA -PAT.7EF-NS .CENTERED. AT LLXL-.-IY 1 ) ANQ IIX^. 

C IY2). THE ANTENNA PATTERN IS STORED AT ANT. DATA M I TH C.25 SPACING 
C AND t AX 6 4 STORAGE ELEMENTS. 3 C TH (IX1.IY1) ANO (1X2. IY2) ARE 

- E- -NE-I-G)-e€HRhGeO-P- Q - iN - T - S - €> f s - X E r~ THfc - C ENTER- OF" THE ANTENN A-F-ATTfcRN.- - FHIS- 
C ROUTINE IS CALLED BY SLBPCLTINE SIGMA. 

C 

-£**************** 1 *********** ****** ******* **************** *************** 

C 


CCMMQN ANT(«4,(4), EM A T I X ( 49 .4 S > .GM A T I X ( 49 ) . P SF ( 6 4 ,S 4 ) 


X l=FLOAT( IX 1 ) 

Y l=FLQ AT ( I Y 1 ) 

X - g»FLOATM » £- ) 

Y2=FLOAT( IY2) 

C 

ETTTEsO.-D 

/■ 

OCL L fl- ~fe=i L..&4 

/* 

Y* (FLO AT ( J >-22. 0)*0.25 

Y 1 M Y -- (Y 1 - Y) 

IF( Y 1MY.LT.-7.7S. OR. Y1MY.GT. S. 0 ) GC TC 10 
Y2MY=-(Y2-Y) 

: I FTY2H Y .CT.-TTTSVOR . VXM YTGTVE.^ J~ GC^TC TO 

J1=32+INT( Y1MY/0.25) 

J2=32+INT< Y SM V>Q.2«) 

C 

dc 11 1 = 1 . e 4 

-c- • 

X* (32 .—FLO A T ( I ) ) *0 .25 
X IMXs- l x 1— x 1 

PPCx ! w X ^CTT- E rOTO H-.~XTMX.O 1 .A./S j C O~~TC~ IT 

X2MX=-(X2-X) 

LELX2-MX.LT £ . Q.QR. X2M X.GT.7 .7S ) . GO TC— LJ _ 

I 1=3 2- INK XJ MX/0.21) 

1 2=3 2— IN T ( X2MXv 0.25) 

—e 

c 

PPOO*ANT( I 1 • J 1 )*ANT( I 2 . J2 ) 

— Ei nZ=EI riSTPROD 

11 CONTINUE 

... -LQ CaNXINUE 

c 

cl 1 1 2=E III Z * C . C 625 


RETURN 

END 


noonn 
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c 

SU8R6U-T- INE - GA MMA-* I S IZ-E-) — 

C 

C***** ****** ************************** 4 ****4*********** *************** 


TO COMPUTE THE GAMMA MATRIX, GAMMA REPRESENTS THE A NT ENN A 
. PATTERN AS. 3URED..ey A JflINDOto .-EUhCTXCJs. _ A .CELT A JEUNCT..ICN IS USED . 

AS ThE AINDUW FUNCTION AND ITS SHIFTING FFCFERTY MAKES GAMMA EQUALS 
TO THE ANTENNA PATTERN ITSELF AT (IX2.IY2). THE GENERATED GAMMA 

€ » AT B I X-- t-S — A — (-1 >— X-t-l-5-I^E E >— €-£-fe-UMN~ M*T«-I X-* 

C 

c* ****** ****** ******************************** ************************ 

COMMON ANT( t4,64 ) .EMA T IX< 49 . 49 ) ,GWATIX<49) .PSF(64 .64) 


DIMENSION IPOINK 49) . JPOINTt 49) 

MSIZE - IS IZE* I SIZE 

C 

ISTX= ISIZE/2 

_ISTY=- ISIZE/2 _ _ _ 

C 

11 DO 10 1=1 .IS IZE 

- IRC iN-T L l L)=-I ST X-T 

1 1=1 I + l 
10 CONTINUE 

jfc- m- x .Lg .MS I ZZ7- CO- TO — M 

C 

_ 11=0 

CO 12 J= L ► I S IZE 
J J= I S I Y-f J— 1 

. OQ-LO—K^ l-I~S IZE. 

1 1= I I + l 

_ JPOINT ( I I ) = JJ 

12 CONT INUE 

C _ 

~ DO ~ 14 T=T » MSIZE 

1X2 = 22 - I NT (FLOAT! I FO I NT{ 2)}/0«25) 

IY-2 -^- 2^-*- INT1FLGA.TLUFCIN-JLI-12-* 0*2S-)- 

GMATIX(I) = ANT( IA2*IY2) 

14 CONTINUE 

RETURN 

END 


i 

i 

i 

-t 


— * 
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c 

SUBROUT I NE COEF ( IS IZE .SN TSR ) _ „ 

C***** ******************** ******** 4*4 ******* ********** *4*44** j* ******* 

C 

C TO iuLVE THE LINEAR EQU"riCN : EVATIX'4 CCEF a G VAT IX. 

C COcF IS THE ESTIMATION COEFFICIENTS FOR GENERATING BRIGHTNESS FE HP- 

C ...ERATURE5- IN TfcE A PC. F-R CCESS* EVATI* ANC GMATIX. ARE CCVRUTcO BY 

C SUBROUTINE SIGMA ANO SUBROUTINE GAVVA. SNTSR IS THE SCUARc OF THE 
C NOISE TO SIGNAL RATIO.* THE RESULTING COEFFICIENTS ARE STORED IN 
C- -GMATIX-REPL ACtNG -THE "GAMMA. — 

************* ****** ************** ***** ******* ************************ 

COMMON ANT(64.e4), EMATIX(49 »49> • G MAT I X (49 ) . P SF (64 .6 4) 

- O-lMEN^iON ..4ICAREAC494 

C 

MSIZE = I SIZE* I S IZE 
DO 10 I* I. MSIZE 

EM AT I X ( I * I ) a EMATIX(I.I) + SNTSR 

10 ~~CCNT TnLE * — ' ‘ 

C 

M '= 1 ' " 

IA a 49 

IOGT a * 

C 

CALL LEQT1FC EMA7IX.M.N .IA .G MATI X . IOGT ,»KA REA .IER ) 

V* 

3UM=0 . 0 

CO — LI- -I? 1 .MSIZE 

SUM=SUM + GMAT IX{ I > 

11 CONTINUE 

S MALL * ( 1 »0-3U M )FM S*ZS 

DC 12 I=t. MSIZE 

GMAT IX ( I ) = GMAT IX( I l+SMALL 

12 ~ CONTINUE 
C 

. RETURN. . 

END 
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— -tHHaRDto TE NE- A NT-P^P— <~ I SIZE) 


C 

C************* ***************** **44**4 4 *** 4*4****-* ** «* * ** * * ** ********* 

C TU COMPUTE THE POINT SPREAD FUNCTION ASSOCIATED WITH A SPECIFIC 

- C._ L^INEAR COMO IN AXIILN — CF_.7H£ ANTENNA -PC LM—SFREAQ FUNCTION ANT _AI X. 

C WHERE PS FIX) = SUMMATION (M 4 ANT(>)>. A IS THE CORRECTION COEFFI- 

C C I ENT MATRIX STOREC IN GMATIX AND A NT IS THE ANTENNA PATTr.fiN. IN 

— e-^-TH 6 R o Q Rgfr, THE P OI N T £ PRE AO — P^N CT I O N » -P-S-F-r— I S - 6B T -* ! NE-C - 6V- APRL-f 

C INC A ISIZFXISIZE LOCAL OPERATOR TC EACH AND EVERT PIXEL OF THE 
C ANTENNA PATTERN. THE LCCAL OPERATOR IS A LINEAR C C PE I N AT ION UPER- 

C’ ATCR - ¥TTF UTVEN CCEFFrciENTS~STXPnT'AT"CFAT IX. THE“RESULTTNG PSF 

C IS STORED ON OISK PSF. DATA (I/O UNIT A 10.) 

C. . _ . ... 

C* ************ ******************** 4** 4 4 *************** ***************** 

C 


- CO MMON -4 


0 I MENS IL N A (49) 


rP^SF-(6 4^» ) 


' TD IFF= ISTZEV J 

ID IFF* ID IFF* 4 

— MSIZE=1SJZEAISIZE. 

C 


DO IQ J= 1 • 64 

„ l KIA nT - 1 T rv f »g 

nt w r^'' , ’T3r * |*r jri 


CO II 1= l * 6 4 
—rsTxRT!=r-rarFF- 
JNDX=JSTART-4 


_K= l 

c 


DO 12 JJ= : . I SIZE 

J N OX » J N O -X*>4- 

IN DX= I ST ART- 4 
C 

• do 12 u = r.mzE 

INOX* INOX+4 

iFl JND X »LE.j« Q .QR . -LNDX.Ja.1 «.&AJ GJO TC LA 

I F ( IN OX .LE . 0 .OP .INDX.GT.64) GC TO 14 
A(K)s*ANT( INOX.JNOX) 

GG-F0 — ^5 — 

14 A(<)=0.0 

15 K=K + 1 


15 CONTINUE 
.12 CO NT INLE_. 


SUM=Q .0 

nn i f, i*— i ucm c 

SUM= SUM+GM AT I X ( M ) *A ( M ) 
16 CONTINUE 

PSF( T,J7»SUM 

C 


11 CONTINUE 

10- CONTINUE 

C 

CO 21 J= 1.64 

WRIT* t t 3 -tPIFtTttT. I~»T 

20 FORMAT (64E 12 .6 ) 

21 _ CONTINUE. 

RETURN 
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SUE RCXTTTNE 0FT2CCM .N »FXY .FU'VRL .FUVIM ) 

C 

C *_******** ************************** * ********************************** 

/■ 

C TWO 0 IMENS IONAL C1SCREIE FCURIER TRANSPCFF. 

-C — FXY— RE PRES EN TS ■ -AM-M-*- N- AFFAY- O F- I N-P UTS- HBE-AL > Cl* A- GP-IC-fiP- *4— A — N 

C POINTS. EACH NUMBER IN THE ARRAV REPRESENTING THE VALUE CF THE 
C FUNCTION AT THE CORRESPOND ING SAMPLING POINT IN Th* X-Y SPACE. 

C FUVRL »' AN" W~ X"N ‘OUTPUT MArTRTrr TS THE BEAT. COMPONENT CF THE DISCRETE 
C FOURIER TRANSFORM OF FXY. FUVIM CONTAINS THE IMAGINARY COMPONENTS 

C OF __TJH.E_TRAN..SFORM.« 

C 

C****« ******** ********** **** ********* A ****** 4 ♦♦♦**♦*** ***************** 

0 I MENS ION FXYO .M .FUVRLO* »N) .FLVIMN .N) 

CATA PJ/^M i ; ?2jC5/ _ 

MN=M*N 

1ST AR T=r_M/! 2 _ - 

JSTART=-N/2 

C 

O*— tOO-r±-t»+* — * 

C 

V=FLO AT! J STAR T +J } 

VN*V /FLOAT (N T 

C 

_ _.DO .ZO(L I=L.M - ... 

C 

U=FLQ AT! ISTART+I ) 

UM=U /-FL&A-T4 M-) 

SUMRL= 0. 0 

SUMIM*C.O _ 

CO 10 JJ^l.N 


Y= FLOAT! JJ-l ) 

*UN=.v*VN.__ 

C 

00 20 1 1 * l • M 

* »F L OATH 1 t—f!r — - 

TH£TA«2.*P I*! X*LM+YVN) 

SUMRL* SUMRL+FXY! I1.JJ)*CCS! THE TAJ 
SUM I Mi SUM I M— F X > (II. J J )'*S I N ! THETA ) 

C 

-- 20— CONT-UUJE — - 

10 CONTINUE 
C 

FU *tt L -fr J N vU -N»SOMfHC/Mfr 

FUVIM! I. J ) = SUMIM/MN 

C_ 

' 200 CONTINUE 
10 0 CONTINUE 

~ RETURN ~ 

END 
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— - St>£ROU-TTN£- OFT 2-XYi f*XY) -- 

► ** ********************** 4 *****4** « 4 4 4************* * t*M**U**n««* 

( 1 ) TO COMPUTE 2-0 DISCRETE FOURIER TRAHSFCFH CF FX Y (64,64.) ALONG 
Tfc£. U- AXIS. AND THE V-AXI S CENTERED AT CQ.d}* 

(2) TO COMPUTE THEIR MAGNITUDES ANC CB VALLES. 

(3) TO PRINT OUT THE RESULTS. 


*♦**♦♦’*•**♦♦♦* *4«*4> ♦*♦♦♦*♦ ♦♦♦♦♦♦♦** ♦♦♦♦Y4444444 ******** ****•♦♦**♦**♦♦ 

~ 'C7METT3T0N FXYf<4.'f 4 ) .FLREALC t4r.FXIMAG(64) ,FVREAL{64) ,FVIMAG<64>. 
♦ FUMAGNt £4 > i,FVMAGNC £4) ,FUOE C64) ,RVD8<64> 

DATA P 1/3.14 l£<2d« y 

<<<<<<<<<< CFT ALONG <U»V) = (U.O) D 5 3 1 0 3 0 5 J 3 DD 


DO 200 1*32, €4 
U=FLO AT( — 3 2+ 1 i 

o«=u/e4. - — - ’ 

SUMRL* O. 0 

. sumim* a^a 

c 

CO 10 JJ*1.«4 

DU 20 1 1=* l » £ 4 
XsFLOATt Il-li 

THETA=r.O*PT4rx*UM J 

SUMRL=SUMP.L+FX>( II ,JJ MCCS< THETA) 

SUMIMF SUM IMr*_F_XY LI I.*. JLLL# SlJVLT HE.TAJ _ 

C 

20 CONTINUE 

HO — •e&HTIHfec 

C 

FUREAL C I )*SUMRL/ 16. 

FUTM AGi r j«STJMTNr7IEi 

C 

. £00 _ CQNlimiE. 

c 

C<<<<<<<<<<<< OFT ALONG (U,V0 * (C.V3 3 3 D 3 3 3 3 > >>» 


C 

C 


DO 100 J=32,£4 
V=FLQAT(-32+J ) 

- vn^v/c-4. 

SUMRL* C • 0 

SUKIM*a^jQ 


DO 30 00*1,64 


YVN* Y*VN 


D O "4. O ' 
THETA* 


£.0*PI*< YVN) 


GK>u*i 

OF POOR QUAUTY 


SUMRL=£U.vtKL+FX>( I I , JJ ) *CCS(T)-ETA) 
^iUM-L'4=3JM.lM^FXYC 11*04 H i.LN(. TfcE-TA-1 

40 CONTINUE 

££, *jT | NL C 

c 

FVfii-u (J >=SUMRL/ 16* 

FVTMAGa >=SUMr»/'167" 

c 

_100 _.CCNTlfUJE . . 

/• 

C<<<<<<<<<<<< COMP LIE MAGNITUDE ANO Oe >>> 3 > XD ;»>>>> 


00 JOG 1=32. £4 

FUMAGNtl )— SQRT ( FUT3AL ( I ) *F UREAL ( I ) +r LI MA G ( I ) *FUI VAGI I )) 
FVMA JN ( r)=STQPT(F YkEAITC rT*FVSEAU fl )>F VTVTG { I) *FVI A AG (I )) 
FUC8I I 1=20 .0*ALUG IC(FUVAGN(I ) ) 

£.V£.i (U= 20.a*ALaGL!l(F.Yl£AGNX2 11 _ 

300 CONTINUE 
C 

C <«■ — G fc-T— »»>»*->»-> 

C 

*r i rg( e. 43 ) 

AT FORMAT ( /~) ■- * — 


l BMAftdl ~l 1 .26.22 - VOU=SCROCO. 0SN=VCK JE. F CLA AC .FCPT 


*R ITE ( £ .44) 

44 FORMATS ' ’ 0=A3<TS TPFAC7 — 

«RITE ( £ • 45 ) IFUREAL! I ) .1=32. £4) 

45 EilRMAT L7L !£_*.£ 1 0 «E.l 1 

-KITE ( C. 56) 

56 FORMAT!/.* TFE O-AXIS (IMAGINARY) :*./) 

»R1-T5-(-6tA6 ) (F UtEMt MH 1 > 

»R ITS ( C . 43 ) 

4d FORM ATI/. ' TFE UH4MS (MAGNITUDE) :*./> 

FdR ITE ( 6.' _ 457‘ CFUMAGNTTT »T"-32V€ - 47 

•RITE(£. 49) 

49. FORM AT! /^*_ TFE_Lr-.AiaS_L_.aE_l 1*-^)- 

*R-iTE(£.45) (FUDE( I ).I=32.£4) 

W -R frT-E t £- . -46-) 

46 FORMAT!/.' TFE V-AXIS (REAL) : * ./) 

mR ITE ( 6. 45 ) ( F VREAL ( I ) . I =32 . £ 4 ) 

~ ~ T AKlTEX£V47r~ " 

47 FORMAT!/.* TFE Y-AXIS ( I VAGI f> Afi Y ) : * ./) 

. JRITEU..451 . CFVIMAG! L)^I=3EuL4J .. ... 

<R ITgC £. SO ) 

50 FORMAT!/. • TFE V-AXIS (MAGNITUDE) : * ./ ) 

— Kfi 1 r c-!6-.-45-)- — fFYMAG N ( I hrl »32 r -*-*-) — ■ — — 

*R ITE! 6. Si ) 

51 FufiMAT!/.* TFE V-AXIS ( C3 ) : * ./ ) 

ARITc( ei 457 CFVCE ( 1 1 • 1 = 22 .64 / “ _ ~ “ 

_ .RETURN - . 

END 
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3SLGhA.^ NftMt : PL C T FCLP.IOF TPaNSFCFF { PLC TFT) 

TJ COMPUTE OF 7 OF AM. Per (47 JX3, 5*5, 7X7). ALL AT J=u. ANO V 

J IM_N£ ICN Ai>, TC t 4 . c 4 ) .P «P { 64 . c'4 > 

Cl. -ad J-- 1 i 6 4 

Kt Ac < i i . d 1 ) (AnT(I,J),Is1,64) 

- ^icAl ( ia. 61 ) (F^Sf-( I ,J ) .1=1 ,t4 ) 

d l r L A,-! AT(t*teit.t) 
cdU s_Li'-»T T ►. Li 

if » i £ i ( i dj ) 

FjAilATl/. 1 jF T CF ANTENNA **♦♦*#****♦*»,>') 

C <>Li_ 3F~2 aV( ANT ) 

IT fc ( t . S 1 Y - - — — 

■;i FCS;iAT(/. • *4*$$4->*** OF” CF FCINT eFAcAC FUNCTION ********* * », / 

L AL*_ ;. C 72XY!P;F ) 

is TOP 
ENO 
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■SEX' 
00010 
00020 
00030 
00040 
00030 
000*0 
00070 
00080 
00090 
00100 
00110 
00120 
00130 
001*0 
00150 
001*0 
00170 
00180 
00190 
00200 
00210 
00220 
00230 
002*0 
00230 
00260 
00270 
00280 
00290 
00300 
00310 
00320 
00330 
003*0 
00330 
0 1)360 
00370 
00380 
00390 
00*00 
00*1 0 
00*20 
00*30 
00**0 
00*30 
00*60 
00*70 
00*80 
00*90 
00300 
00310 
00320 
00330 
003*0 
00330 
00360 
003T0 
00380 
00590 
00600 
00610 
00620 
00630 
006*0 
00630 
00660 


ORIGINAL PAGE IS 
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C PROGRAM NAME * FSCENE 

C TQ COMPUTE the FOURIER TFANtFORM OF AN IMAGE USING SUBROUTINE 
C DFT2D. THE OPIGINAL IMAGE IT STORED AT SCENE. DATA WITH LUN-1*. 

C THE RESULTING PEAL AND IMAGINARY PARTS ARE STORED AT F'CENEP. DATA 
C AND r SCEhEI . DATA. WITH 1? «- 16 At THE LUN' S- 


INTE6ER*2 IDUM<6*.6*> 

~ T MEMS I ON FXY < 6*« 6*> • BUFFPL <6*» «BUFFIM<6*> 


N»6* 

LUN-1* 

LUNRL-13 
LUN I H— 16 
REWIND LUN 
C 

DO 10 J-l,N 

READ < LUW> < I3UM< I » J> • I —1 ,N) 

DO 11 I— 1 »M 

FXY ( 1 1 J> » I DUM < I » J) 

11 CONTINUE 
10 CONT INLiE 
C 

CALL DFT2D ' H* N* FXY'* LUNRL, LUN IN) 

C 

REWIND LUNRL 
REWIND LUNIN 

C 

C< - • <•. READ and PRINT OUT ^ , 


WRITE ’6* 1 13> 

113 FORMAT < ' PRINT OUT INAGE cYES-l) ' • 
READ <3, 11*) IPRINT 
11* FORNAT »' 1 1 > 

IF c I PRINT. ME. 1) GQ TO 100 

A-l. 0 
IE-0 


DO 13 INDX— 1 »* 

IS-IC-l 

IE-IS+13 

DO 1* J— 1 » 6* 

READ <LUNRL> (BUFFPL < I) * 1-1 »6*> 

READ <.LUNIM> <BUFFIM« I) . 1-1 , 6*) 

DO 13 I-IS-IE 

BUPFRL C I) -<SQPT <BUFFRI_(I> *BUFFPL * I) ♦BUFF IN f I > *BUFFIN^ !:■ ) > *A 
13 CONTINUE 

WRITE C6» 112) <BUFFRLcI),I«IS,IE) 

112 FORNAT C16F3.2, 1 • 

l* CONTINUE 

REWIND LUNRL 
REWIND LUNIN 
WRITE C6» 123) 

123 FORNAT (.Sts'* 

13 CONTINUE 

100 STOP 
END 
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